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Abstract. We consider quantum random walks (QRW) on the 
integers, a subject that has been considered in the last few years 
in the framework of quantum computation. 

We show how the theory of CMV matrices gives a natural tool 
to study these processes and to give results that are analogous 
to those that Karlin and McGregor developed to study (classical) 
birth-and-death processes using orthogonal polynomials on the real 
line. 

In perfect analogy with the classical case the study of QRWs on 
the set of non-negative integers can be handled using scalar valued 
(Laurent) polynomials and a scalar valued measure on the circle. 
In the case of classical or quantum random walks on the integers 
one needs to allow for matrix valued versions of these notions. 

We show how our tools yield results in the well known case of 
the Hadamard walk, but we go beyond this translation invariant 
model to analyze examples that are hard to analyze using other 
methods. More precisely we consider QRWs on the set of non- 
negative integers. The analysis of these cases leads to phenomena 
that are absent in the case of QRWs on the integers even if one 
restricts oneself to a constant coin. This is illustrated here by 
studying recurrence properties of the walk, but the same method 
can be used for other purposes. 

The presentation here aims at being sclfcontained, but we re- 
frain from trying to give an introduction to quantum random walks, 
a subject well surveyed in the literature we quote. For two excellent 
reviews, see [HUH]. See also the recent notes [19] . 
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1. Introduction and contents of the paper 

We start with a brief look at classical random walks. 

Consider, for simplicity, a discrete time random walk in a denu- 
merable state space which we take to be the non-negative integers 
i = 0, 1,2 

The state of the system at time n is given by a row vector -K^ n . This 
is also called the probability distribution at time n. The i th component 
is interpreted as the probability that a particle can be found at site % 
at time n. 

These non-negative quantities 7Tj <n are assumed to add up to one, 
for any n, when summed over i in the set of non-negative integers. 
This is not the only way to describe these simple processes, but this 
facilitates the transition to the quantum case. 

The evolution of the system is given by some transition probability 
matrix P = (i\j). This means that for each state i we have a collection 
of "transition probabilities" Pij with the condition that these non- 
negative numbers add up to one when summed over the index j. P it j 
is interpreted as the probability of a transition form site % to site j in 
one unit of time. 

The Markov nature of our process is given by the fact that at time 
t = 1 the new probability distribution gives weight 

oo 
j=0 

to the i th site. In other words the state at time t — 1 is the vector 
obtained by multiplying the probability distribution at time 0, namely 
TTjo, by P. The Markov property says that the process starts from 
scratch at any integer time, so that the probability distribution (i.e., 
the state of the system) at time t = n is obtained by taking the product 
of the row vector 71^0 and the matrix P n . 

Random walks with a stationary (i.e., time invariant) transition 
mechanism as above have been studied extensively. In some cases it 
is convenient to exploit some of the mathematics associated with the 
matrix P. We will not go into details here but it turns out that this 
is the case when the matrix P is either symmetric to begin with, i.e., 
Pi j = Pj A or it is what is called symmetrizable, namely there exists a 
vector 7Tj such that 

This notion is on the one hand related to the issue of reversibility 
and on the other it allows one to introduce a certain inner product in 
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L 2 (Z > 0) that makes P selfadjoint and thus one has recourse to a nice 
and simple spectral theory. 

A new item comes in now: the idea of a "local" transition. More 
specifically we consider a special kind of random walks, namely those 
that only allow for "nearest neighbour" transitions, i.e., Pjj = if 
the indices i,j differ by more than one unit. In this case the matrix 
P is tridiagonal and this is the best of all worlds since under very 
simple conditions any such matrix is symmetrizable into a tridiagonal 
one. Here we see a physically important issue such as nearest neigh- 
bour interactions going along with a nice mathematical fact such as 
symmetrizability. 

The consequence of the combination of these two features is that 
one has a powerful and natural tool to study random walks with a 
tridiagonal or Jacobi transition matrix, namely the classical theory of 
orthogonal polynomials on the real line, a subject that goes back to 
the beginning of the 20th century. 

This idea of making explicit use of the spectral theory for selfadjoint 
operators in Hilbert space to study an important class of random walks 
on Z > probably appears first in [IB] . The authors point out that 
similar ideas had been used earlier in the case of diffusion processes by 
W. Feller, see [ID], and H. P. McKean, Jr., see [23]. One could add to 
this list of precursors of this fruitful line of work other papers, such as 
[22] . One can also handle the case of continuous time, but we do not 
pursue this here. 

The basic idea in [16] is simple: if one is dealing with these Markov 
chains that only allow for a "local" transition with Pjj = if i, j differ 
by more than one unit (also known as birth-and-death processes) then 
the orthogonal polynomials that can be built out of the corresponding 
Jacobi matrix and the underlying orthogonality measure give all the 
ingredients of the spectral resolution of the matrix. In particular they 
allow for a simple way to compute its power P n . 

All of this can be adapted to the case when our state space is not 
the set of non-negative integers but the set of all integers. The fact 
that now we have two different ways of going to infinity calls for the 
use of matrix valued orthogonal polynomials, a notion introduced by 
M. G. Krein [201 EI] around 1950. In this case we need polynomials 
with values in the set of square matrices of size two, whereas up to now 
we were dealing with scalar valued polynomials. 

In fact a larger class of random walks, known under the name of 
quasi-birth-and-death processes, can also be analyzed by using these 
matrix valued orthogonal polynomials on the real line. In many cases 
the orthogonality measure is not really a matrix valued measure, but 
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an appropriate linear functional. For two papers doing just this see 

002]. 

With this background we can now move closer to the subject of our 
paper, namely the study of a different sort of processes that go under 
the name of "quantum random walks" (QRW). 

The basic idea is this: a system has a (denumerable) set {\i)}iei of 
measurable states called "pure states". The state \ty n ) of the system at 
time n is a (complex) superposition |\l/ n ) = £) ier V^nl*) °f P ure states, 
so it is described by a "wave function" fy >n with % running over the 
pure states. Therefore, a state can be identified with its wave function. 

The complex number ipi >n is no longer a probability but a probability 
amplitude, so that the actual probability to be in the state \i) at time 
n is l^nl 2 - The total probability of finding the system at time n must 
be 1, thus the wave function must satisfy the condition 

(1) £IW 2 = 1 - 

iei 

In other words, the wave function lives in the unit ball of L 2 (X), which 
is equivalent to consider the pure states as an orthonormal basis of a 
Hilbert space whose unit vectors are the states of the system. 

According to the conservation of the probability, in the time invari- 
ant case the time evolution of the system is characterized by a single 
unitary operator il on the Hilbert state space: \^ n ) = H n \^ ). If 
U = (Uij)ij £ j is the unitary matrix such that H\i) = X^ez^»>jb')> ^ e 
evolution of the wave function is governed by 

(2) = = (V'i 1 n)ieX- 

That is, the evolution operator in the wave function representation is 
a unitary operator U on L 2 (I). 

When X = Z, it is proved in [23] that all the local QRWs (U itj = if 
\z — j\ > r for some r) which are translation invariant (Uij = f/j+ij+i) 
are, up to a phase, integer powers of a simple translation, i.e., il = 
e *T»,T\j) = \j + l). 

If we do not ask for translation invariance but allow only nearest 
neighbour transitions (C/ 4 • = if \z — j\ > 1) then the QRW splits 
into independent QRWs with no more than two pure states or, up 
to a change of phases of the pure states, it is a simple right or left 
translation, i.e., = e l0J \j ± 1). The situation is even worse in the 
case I = Z > 0, where the restriction to nearest neighbour transitions 
always forces the splitting. These results, less known in quantum com- 
putation, are direct consequences of [H Lemma 3.1] and [U Theorem 
3.9]. 
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A way to generate non-trivial QRWs on the integers or the non- 
negative integers, allowing at the same time only nearest neighbour 
transitions, is to include extra degrees of freedom. Following the anal- 
ogy with a quantum physical system, the simplest way to do this is to 
consider that every site i 6 I has an internal degree of freedom of spin 
1/2 type, i.e, taking values s G S = {j, J.}. The space state is a tensor 
product spanned by {\i} <g> \i) <g) ||)}j 6 z- 

We will be considering mostly the cases when X is either the set of 
non-negative integers or the set of all integers. The evolution of the 
wave function ipi tS>n satisfies ([[J when summing in i G X and s G S, 
so its evolution is given similarly to (j2J) by a unitary operator U on 
L 2 (lxS). 

Although the choice X = Z > would be as natural as X = Z, 
the first one has not received too much attention in the quantum case, 
maybe due to the difficulties to work with a non translation invariant 
system. However, just as in the classical case, QRWs on the non- 
negative integers are more natural for an orthogonal polynomial ap- 
proach. Indeed, from this point of view, QRWs on the non-negative 
integers will be the cornerstone for the analysis of QRWs on the inte- 
gers. 

We are finally ready to state the purpose of this paper. We will 
show that in the case of a large class of quantum random walks on 
the integers there is a natural tool that takes the place of the ma- 
trix valued orthogonal polynomials on the real line briefly alluded to 
above, namely the theory of Laurent 2x2 matrix valued orthogonal 
polynomials associated with a certain kind of unitary matrices, namely 
CMV matrices. This is a ready made tool that combines the necessary 
features for a quantum mechanical description of the phenomenon of 
nearest neighbours transitions: unitarity and a block tridiagonal shape. 

In the next section we first review the standard tools to deal with a 
simple random walk on the integers: the first two are very well known 
and have been used in the case of quatum random walks. The third 
method, involving matrix valued polynomials on the real line is less 
well known as a useful tool for the study of classical random walks on 
the integers. To the best of our knowledge neither scalar nor matrix 
valued orthogonal polynomials have so far been used in the case of 
QRWs and they constitute the main novel point of this paper. 

The contents of the paper are organized as follows: 

Section 2 reviews different approaches to the study of classical ran- 
dom walks, with special emphasis in the methods using orthogonal 
polynomials. This will establish a benchmark for the development of 
analogous techniques in the quantum case. 
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Section 3 introduces CMV matrices and Szego polynomials. 

Sections 4, 5 and 6 introduce QRWs and consider two extremely 
simple examples to show the workings of our method. 

Section 7 considers the case of a not (necessarily) constant coin. 
We refer to this as dealing with "distinct coins". In principle our 
method can handle this general case, something that the more standard 
methods cannot do. 

Section 8 takes up the simpler case of a constant coin and one sees 
that the case of the integers reduces to the study of two QRWs on the 
non-negative integers. 

Section 9 tackles in detail the case of a QRW with a constant coin 
on the non-negative integers. The study of this case will be used in the 
following section and it is of independent interest since the phases of 
the coin have now a strong influence on the results, and in general the 
orthogonality measure has a discrete mass. 

Section 10 uses the results of the previous two sections and shows 
that when dealing with all the integers there is no mass point in the 
orthogonality measure and that the phases of the coin play no role. 

Section 11 applies the results of the previous sections to compare 
in the case of two specific QRWs their behaviour when considered on 
Z > or Z. 

Section 12 is devoted to getting some large time asymptotics. 

Section 13 takes up the issue of recurrence of a QRW. By using our 
method we see here some marked differences between the classical and 
the quantum case. This section exhibits very clearly the benefits of 
using the notions introduced in this paper to analyze QRWs. 

Section 14 lists some conclusions as well as some open problems 
that can be treated with the methods in this paper. 

2. Classical random walk on the integers, three looks at 

a classical subject 

Consider a random walk on the integers with p, q respectively the 
probabilities of going right or left in one unit of time starting at any 
integer position i. 

Denote by P(a, b, n) the quantity of interest, namely the probability 
of going in n steps from the initial position a to the final position b. 

We describe three ways to study this basic problem. 

a) Path counting: we must have c steps to the right and d steps 
to the left with c + d = n and c — d = b — a. One can thus solve 
for c, d in terms of a — b, n. Each such path has probability p c q d 
and there is total of (™) such paths. This gives P(a,b,n). 
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b) Fourier methods: if for simplicity we take p — q — 1/2, it is 
clear that 

P(a, b, n) = ^{P{a, 6 - 1, n - 1) + P(a, b + 1, n - 1)). 
Introduce the Fourier series 

oo 
b=— oo 

with an inverse given by 

P(a,b,n) = ±- ^ e- m P:(6)de. 
27r Jo 

The difference equation above can be used to see that 

P2{o) = pr\&) cos e = pM(coser 

and therefore 

P(a,6,n) = — / P° (6) (cos 6) n e- lbe d6. 
In the special case when P(i,j, 0) is given by 5y we get 

P(a,b,n) = — f (cos9) n e- i{b - a)e d9. 
2?r Jo 

This integral can be computed explicitly in terms of Bessel 
functions, or just as easily, it can be seen to agree with the 
general expression obtained earlier. For general values of p, q 
one simply replaces cos# by a weighted sum of e* 61 , e~ tB '. 

These two methods have been properly adapted to study 
QRWs on the integers. "Path counting" was used by D. Meyer, 
see [23]. This analysis was pushed further along by using Ja- 
cobi polynomials in (2]. In this same paper as well as in [25] 
one finds expressions obtained by using "Fourier methods" to 
analyze the appropriate recursion relations. The reader can 
find a very sophisticated use of these formulas to derive some 
aymptotic results about the Hadamard QRW in [7]. A recent 
and very precise analysis of asymptotic results is given in [T^j . 
Some interesting use of generating functions to obtain rigorous 
results is made in [3]. 

c) Using matrix valued orthogonal polynomials: since this is not 
too standard and since this is the way in which we analyze 
the quantum case we give a lengthier account of this way of 
tackling this problem. This method goes back, in spirit at least, 



8 



CGMV 



to the original paper of Karlin and McGregor, |16j . where the 
reader will find a complete discussion of the case of a random 
walk on the non-negative integers using scalar valued orthogonal 
polynomials. For a fuller account of the material below see 

[SUM- 

Consider, with M. G. Krein, [20j 12], the set of polynomials of 
the real variable x with matrix coefficients of a fixed size d. All the 
matrices that appear below have this common size d and this is the 
appropriate choice when the state space is the cartesian product of the 
set of non-negative integers with the set 1, 2, 3, . . . , d. 

Given a positive definite matrix valued measure dM(x) with finite 
moments j R x n dM(x), n = 0,1,2,..., consider the skew symmetric 
bilinear form defined for any pair of matrix valued polynomial functions 
P(x) and Q(x) by the numerical matrix 

(P,Q) = (P,Q) M = [ P{x)dM{x)Q{x)\ 

where Q(x)^ denotes the conjugate transpose of Q(x). 

By the usual Gram-Schmidt construction this leads to the exis- 
tence of a sequence of matrix valued orthogonal polynomials P n {x) = 
K n , n x n + K n:n -\X n ~ Y + • • • with non-singular leading coefficient K n ^ n . 
We make no special assumption on K n ^ n . 

Given an orthogonal sequence {P n (x)}n>o of matrix valued orthog- 
onal polynomials one gets by the usual argument a three term recursion 
relation 

(3) xP n (x) = AnPn^x) + B n P n (x) + C n P n+1 (x), 

where A n , B n and C n are matrices and the last one is non-singular. If 
we had insisted on orthonormal polynomials then we would get some 
relations among these coefficients. 

It is convenient to introduce the block tridiagonal matrix P 

/B Co \ 
P = A B x C x 

If Pi j denotes the i, j block of P we can generate a sequence of d x d 
matrix valued polynomials Q i (x) by imposing the three term recursion 
given above. By using the notation of the scalar case, we would have 

PQ{x) = xQ{x) 

where the entries of the column vector Q(x) are now d x d matrices. 



MATRIX VALUED POLYNOMIALS AND QUANTUM RANDOM WALKS 9 



Proceeding as in the scalar case, this relation can be iterated to 
give 

P n Q{x) = x n Q(x) 

and if we assume the existence of a positive definite matrix valued 
measure dM(x) as in Krein's theory, with the property 

(Qj,Qj)6 id = / Q i (x)dM(x)Q j (x)\ 
Jr 

it is then clear that one can get an expression for the i, j entry of the 
block matrix P n that would look exactly as in the scalar case, namely 

(P n ) isj (Q j ,Q j )= [ x n Q l (x)dM(x)Q J (xy. 
Jr 

The expression above, allowing one to compute the entries of P n is 
usually called the Karlin-McGregor formula, see [T6j . 

It may be worth noticing that the integrals above are different from 
the ones that appear when one uses the Fourier method. The same will 
be true in the quantum case. 

Just as in the scalar case, the expression above becomes useful when 
we can get our hands on the matrix valued polynomials Qi(x) and 
the orthogonality measure dM(x). When this is the case this formula 
allows us to compute the transition probabilities between any pair of 
states i < j in any number of steps by using only the top j rows 
of the matrix P. The "time dependence" has now been isolated to the 
term x n and the study of its behaviour for large values of n can be 
handled with traditional methods. The same remark will apply in the 
quantum case and this will be used in Section [T2J 

We are now ready to tackle the example of random walk on the 
integers, when the probabilities of going right or left are p and q re- 
spectively. This is the most general translation invariant random walk 
on the integers with nearest neighbours transitions. We present it here 
to compare this analysis with its quantum counterpart, i.e., the quan- 
tum random walks with a constant coin on the integers analyzed in 
Sections M and [TU1 

If we "fold" the integers by relabelling the natural sequence 

•■■-3,-2,-1,0,1,2,3... 

in the fashion 

...5,3,1,0,2,4,6... 

then the transition probability matrix goes from being a scalar tridi- 
agonal doubly infinite one with p in the 2,2 + 1 diagonal and q in the 
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i + 1, i diagonal to the following semi-infinite block tridiagonal matrix 
(with 2x2 blocks) 

\ 



p + q = 1. 



/o 


Q 


V 

















V 








Q 














q 











p 














P 











q 














Q 











p 














p 











q 



V 



For this first example the appropriate matrix measure is already 
found in the original paper by S. Karlin and J. McGregor. It is given 
by 



dM{x) 



1 x/2q\ 

^A P q - x 2 \ x / 2 q p/q ) 



dx, 



\x\ < y^pq. 



One does not find in the paper mentioned above the corresponding 
matrix valued orthogonal polynomials or the block tridiagonal matrix, 
but they can be easily given, see [12] and below. 

In this example we have 



B 



q 
p 



Bi 








k > 1, 



Ai 



q 
p 



k > 1, 



p 
q 



k > 0. 



The orthogonal polynomials given by 

A k P k -i(x) + B k P k (x) + C k P k+1 (x) 
P^(x) = 0, P (x) = 1, 



xP k (x) 



where and 1 are the null and identity matrix, can be easily expressed 
in terms of Chebyshev polynomials. 

Let us denote by U n (x) the Chebyshev polynomials of the second 
kind, which satisfy 

(4) U n+1 {x) + U n ^{x) = 2xU n {x), U-!{x) = U {x) = 1. 
The relation with the Chebyshev polynomials U n (x), is given by 



Pk(x) 



(q/p) 2 o 

(p/q) 1 - 



lUdx*) - 



^ (q/p)' 
(p/q)^ 



U k -i{x*) 



where x* = x/2^/pq. This expression can be compared to the one in 
Section [TUJ 
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These ideas have been used recently to analyze more elaborate ex- 
amples of classical random walks. The reader may consult [HI d21 E321 

03 EES]. 



3. SZEGO POLYNOMIALS AND CMV MATRICES 

All the information about a QRW is encoded in the unitary operator 
U governing the evolution of the system. Therefore, it is not strange 
that the theory of canonical matrix representations of unitary operators 
on Hilbert spaces should play an important role in the study of QRWs. 
Surprisingly, such a theory has been developed only recently (see [51 
[32]) giving rise to the so called CMV matrices, related to the Szego 
polynomials. 

Due to their relevance for the rest of the paper, we will summarize 
in this section the main facts about CMV matrices. 

The basic idea is that, as a consequence of the spectral theorem, 
any unitary operator is unitarily equivalent to a direct sum of unitary 
multiplication operators, i.e., operators of the type 

(5) U»: LJ(T)^LJ(T) 

/(*) — > zf(z) 

H being a probability measure on the unit circle T = {z e C : \z\ = 1} , 
and L£(T) the Hilbert space of /^-square-integrable functions with inner 
product 

(f,g)= fW)g{z)d l x{z). 

Thus, it is enough to discuss the canonical representations of unitary 
multiplication operators. Moreover, we can suppose that fi has an 
infinite support, otherwise L 2 (T) is finite-dimensional, so is unitarily 
diagonalizable. 

Since the Laurent polynomials are dense in L^(T), a natural basis 
to obtain a matrix representation of £7 M is given by the Laurent polyno- 
mials (Xj)'jLo obtained from the Gram-Schmidt orthonormalizalization 
oi{l,z,z-\z 2 ,z- 2 ,...}m LJ(T). 
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The matrix C = (xj, zXk)^k=o °f ^ with respect to (Xj)fLo has the 
form 



(6) C = 



/ QtQ 


Po«i 


PoPi 














Po 


— a^Ui 


-a Q pi 



















—a{ct2 




P2PZ 











P1P2 


-aip 2 


— a 2 a 3 


-a 2 p 3 

















p 3 a 4 


—a 3 a 4 


p 4 a 5 


PAPb 











P3PA 


-a 3 p 4 


—a 4 a^ 


-a 4 p 5 


v... 















where pj = ^/l — \oij\ 2 and (oj)jl is a sequence of complex numbers 
such that \oij\ < 1. The coefficients <x,- are known as the Verblunsky 
parameters of the measure p, and establish a bijection between the 
probability measures supported on an infinite set of the unit circle and 
the sequences in the open unit disk. 

Another equally natural basis would be the Laurent polynomials 
(xj)^2_Q obtained from the orthonormalization of { l,z 1 ,z,z 2 \z 2 ,...}. 
They are given by 

x i( z ) = Xj(i/X> 

and, consequently, the matrix of with respect to (xj)°^ is the trans- 
pose C T of C. 

As a consequence, we have the identities 

x(z)C = zx(z), x = (Xo,Xi,X2,---), Xo = l, 
Cx(z) = zx(z), x = (x ,xi,X2, ■ ■ -) T , x — 1, 

which can be viewed as recurrences which determine the orthonormal 
Laurent polynomials. 

The unitary matrices with the form ([6]) or its transpose are called 
CMV matrices, but we will reserve this name for the matrix in ([6]). 
Also, when talking about orthonormal Laurent polynomials, we will 
refer to (xj)^ , which are the ones that we will normally use. 

The canonical representations of the unitaries are the narrowest 
banded representations that can be obtained for all such operators. 
The previous results state that every unitary operator has a matrix 
representation which is a direct sum of CMV matrices, so the canonical 
representations of the unitaries are at least five-diagonal. That they 
are exactly five-diagonal has been proved in [6], where it was shown 
that not every unitary operator admits a four-diagonal representation. 
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The CMV matrices have also a tridiagonal factorization C = CM., 
with two unitary 2 x 2-block diagonal symmetric factors given by 



£ = diag(9 , 62,64, .. . ), 

6j 



A^ = diag(i,e 1 ,e 3 ,...), \pj 



-a 



The Verblunsky parameters have a special meaning in terms of the 
Szego polynomials (y_j)^=0' wm ch come from orthonormalizing {z J }°^L , 
see [SHUT]. 

These polynomials are not so useful as a basis because the polyno- 
mials are not always dense in L 2 (Y). Nevertheless, they are related to 
the orthonormal Laurent polynomials by 

, g * X2j{z) = z- J tp* 2j {z), X2j+i(z) = z- 3 tp 2j+ i{z), 

x 2j (z) = z~ J ip 2j (z), x 2j+1 (z) = 3~ J ~Vjw+i(*)> 



where <p*j(z) = z 3 ipj(l/z). 

The key result is that the Szego polynomials are determined by the 
recurrence relation 

(10) pjtp j+1 (z) = zifj{z) - ajip*(z), <Po = l, 

so the recurrence for the monic orthogonal polynomials (0j)°5L o is 

(11) (j>j+i( z ) = Z( Pj( z ) -aj<j>i(z), 0o = 1, 



which shows that <x,- = — 0j +1 (O). 

For some measures on the unit circle the Szego polynomials, and 
therefore the Verblunsky parameters and the CMV matrix, are known 
explicitly, see [28]. Other cases can be analyzed following different 
methods. 

For instance, given a sequence of Verblunsky parameters, the spec- 
tral analysis of the corresponding CMV matrix can be used to obtain 
information about the orthogonality measure because it is related to the 
spectral measure of the CMV matrix. Indeed, the support of the mea- 
sure coincides with the spectrum of the CMV matrix, the mass points 
zq being the eigenvalues, which are simple and have eigenvectors given 
by x(zq). Bearing in mind (j7|), this means that x(zq) G L 2 (Z > 0) 
exactly when z is a mass point. Actually, p({z }) = 1 / 1 1 or ( ^ ) I j 2 - 
These results will be of interest later on, see Section [121 
Perturbative results are useful to study new examples taking as a 
starting point known ones. The simplest example of this is a rotation 
of the measure by an angle i.e. 

dp(z) — > dp(e~ l,d z). 
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The change of the monic orthogonal polynomials 4>j(z) — > e y 'Vj(e _l ' 3 z) 
shows that the effect of the rotation on the Verblunsky parameters is 

(12) aj -> e-tWaij 

while (J5J) implies that the orthonormal Laurent polynomials transform 

as 

(13) x 2j ^(z) e-^.^e-^^), z 2i (*) -> e«%( e - w z). 

The transformation above will play an important role in later sec- 
tions. 

Another tool for the study of Szego polynomials is the Caratheodory 
function F of the orthogonality measure /x, defined by 

(14) F(z)= [t±±dn(t), \z\<l. 

J j t — z 

F is analytic on the open unit disc with McLaurin series 

oo „ 

(15) F{z) = 1 + 2 Y'/ZjV, n,= I z?dn(z), 

3=1 h 

whose coefficients provide the moments \Xj of the measure \x. 
F can be obtained as 

(16) F(z) = lim z < 1, 

i-oo p*(z) 

where </5j are the Szego polynomials whose Verblunsky parameters are 
given by — aj if the original ones were aj. 

The Caratheodory function is a shortcut that allows one to recover 
the measure from the Szego polynomials. If 

d9 

(17) djj,(e id ) = w(6) \-dfi s (e i9 ), 6 G (-7r,7r], fi s singular, 

the weight w(6) is given by 



(18) w{6) =\imReF{re ie ), 

rfl 

and the support of /a s lies on {e 10 : lim r ji F(re ) = oo}. In particular, 
e l9 ° is a mass point of /i with mass /i({e ie °}) if and only if 

(19) m^F(re*)^0. 

Two natural extensions of CMV matrices are of interest to us: dou- 
bly infinite CMV matrices and block CMV matrices. 
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Given a two-sided sequence (aj)j e z we can define the doubly infinite 
CMV matrix C = CAi, where 

£ = diag(. . . , 6_ 4 , 6_ 2 , @o, ©2, ©4, • • • ), 

M = diag(... J e_ 3 ,e_i J e 1 ,e 3 ,...), 

and Qj, given in (jSJ), acts on the indices j and j + 1. We will use the 
same notations as for semi-infinite CMV matrices. 

C is a five-diagonal doubly infinite unitary matrix with the form 



P- 3 a- 2 -a_ 3 «_2 P-2«-i 
P-3P-2 -a_ 3 p_2 -a_ 2 a_i 













P-iPo 



p_ 2 P-i 

-CK-2P-1 



-a_iao 
-a-iPo 



Po«i 
-«o«i 



PoPi 
-«oPi 



Doubly infinite CMV matrices are also related to the other gener- 
alization of interest: block CMV matrices. They appear in connection 
with matrix valued Szego polynomials, see [8l |28| [29]. 

Given a positive definite dx d matrix valued measure \x on the unit 
circle we can define the right and left "inner products" 



g{z)dp{z) f(z) 



(f,g) 



R 



for /, g with values in the set of d x d matrices. Notice that the symbol 
' which denotes the adjoint matrix includes the conjugation of z. In 
what follows we suppose that j T dp(z) = 1 is the d x d unit matrix. 

Now we can consider right and left matrix valued Szego polynomi- 
als, (vj^j^o an d (vfyjLoi arising from the standard orthonormalization 
of {1z^}jLq with respect to (, ) R and (, ) L respectively. 

Analogously to the scalar case, the matrix valued Szego polynomi- 
als satisfy a recurrence given in terms of a sequence (ctj)'jL of d x d 
matrices such that ||cXj|| < 1. However, this recurrence mixes the right 
and left polynomials in the following way 



t,„-R* 



zipj [z) -a'jif 
z<p?(z 



,~ R 

^0 



where <fij(z) = z^tp^l/zY and pj, pf are the positive definite matrices 



1 



«K) 1/2 , 



Pf 



This mixture between right and left polynomials also appears in the 



connection with the orthonormal Laurent polynomials (x 



jjj=0, {Xj)j=0 
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given by 
(20) 



x 2j (z) = z- j cp^(z), x 2j+1 (z) = z~ j ~ l ^* +1 { 



z 



and related by Xj(z) = x J -(l/z')L (x^fLo comes from the orthonormal- 
ization of {1, lz, lz -1 , lz 2 , lz~ 2 , . . . } with respect to (, ) R , and (a^-)^ 
from the orthonormalization of {1, lz -1 , lz, lz~ 2 , lz 2 , . . . } with re- 
spect to (notice a slight difference with respect to [5], where the 
two kinds of Laurent polynomials discussed are both orthonormal with 
respect to (, ) R ). 

The matrix C which determines the orthonormal Laurent polyno- 
mials through the recurrences 

x(z)C = zx(z), X = (X0.X1.X2,---)> Xo = 1, 
Cx(z) = zx(z), x = (x ,x 1 ,x 2 , . . .) T , x = l, 



is given by C = CAi with 

C = diag(0 o , 2 , ©4, ■ ■■). ( a \ 



M = diag(l,©i,© 3 ,...), \Pj -a 

The unitary matrix C is called a block CMV matrix. Its explicit 
form is 



/ c*l p%at\ pLp[ 

p§ -ct at[ -a pf 

ct 2 pf —ot\(Xi 

p^pf -pfai 


pfp§ -p?a 3 -a A a\ -a 4 p% 



























P2<*3 


P2P3 










-«2P3 








«4Pf 


-oc\a 3 


P4«5 


P4P5 



where the symbol stands for the d x d null matrix. 

Particularly simple is the case of diagonal Verblunsky parameters 
ctj = diag(a], . . . , a^). The corresponding matrix orthonormal poly- 
nomials and orthogonality matrix measure are diagonal too. In other 
words, the recurrences in (f2~Tl) split into d ones associated with scalar 
CMV matrices with Verblunsky parameters 

Finally, the matrix valued Caratheodory function 

(23) F(z)= ft±l d p{t), \z\<\, 

Jj t — z 

and its real part KeF(z) = (F(z) + F(zy) /2 allows one to recover the 
matrix measure p just as in the scalar case. In particular, the matrix 
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moments fij of /x come from the McLaurin series of F, i.e., 
(24) F(z) = 1 + 2 V fi]z j , ^ = [ zH»{z). 

4. Quantum random walks 

We will consider a few one-dimensional quantum random walks with 
pure states \i) Cg> ||) and \i) ® ||), where i runs over the non-negative 
integers or over all the integers, and with a one step transition mecha- 
nism given by a unitary matrix U. Our goal here is to give an explicit 
expression for the entries of the matrix U n for n = 0, 1, 2, . . . describ- 
ing the time evolution of our process. More precisely, we will find a 
Karlin-McGregor (KMcG) formula for such entries. 

The relevance of the CMV matrices, as the canonical representa- 
tions of the unitaries, is clear. However, the reduction of an arbitrary 
infinite unitary matrix to CMV form can be a difficult task. Neverthe- 
less, the QRWs usually considered in the literature include only nearest 
neighbour transitions. In this case it is possible to prove that the CMV 
form is obtained after a simple change of phases in the basis, see [HI 
Theorem 3.2 and Lemma 3.7]. Then, the relation of the CMV matri- 
ces with the Szego polynomials provides the functions and measures 
needed to obtain a KMcG formula. 



5. A VERY ELEMENTARY example 

A simple QRW on the non-negative integers corresponds to the 
unitary operator 

oo oo 

(25) 5 + = ^|z + i)(z|®|T)(TI + ^K-i)(^l!)(!l + |0)(0|®|T)ai, 

i=0 i=l 

which means that the transition mechanism moves spins at positions 
0, 1, 2, . . . according to the following (deterministic) prescription: spins 
up at any location move one step to the right, spins down at 1, 2, 3, . . . 
move to the left, and finally a spin down at location reverses orien- 
tation and stays at location 0. In the next unit of time this spin will 
move to location 1. 

If we choose to order the pure states of our system as follows 



(26) |0)®|T),|0)®|1),|1)®|T),|1)®||), 
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then the transition matrix is 








1 














1 
































1 











1 



































1 











1 












We refer to the pure states, ordered in this fashion, as the zeroth, the 
first, the second state, etc. 

However, there is another way to look at the transition matrix. We 
could think of the functions 1 2 Z j Z y Z 5^ 

, . . . defined on the unit 
circle and form a column vector x(z) with these functions. Applying 
our matrix U + to this vector is the same as multiplying the vector by z, 
which is another way of saying that for each value of z we have found 
a formal eigenvector with eigenvalue z for our matrix. 

If we denote the components of this column vector by Xj(z), j = 
0, 1, 2, . . . , we claim that the probability amplitude of going in n (pos- 
itive or negative) units of time from a pure state j to a pure state k 
(the indices j, k run over all the zeroth, first, second, . . . pure states 
introduced above) is given by the integral 



It is completely elementary to see that this integral is either or 1 and 
that in fact we get the correct expression for the amplitudes. 

The above relation is not an accident. In fact, this example is the 
prototype of everything that follows. The matrix given above is the 
simplest example of CMV matrix, corresponding to null Verblunsky 
parameters. The Laurent polynomials (£j)°L are orthonormal with 
respect to the Lebesgue measure dz/2niz = dQjl-n. Therefore, the 
identity Ux(z) = zx(z) gives 



6. Another example 

The natural extension of the previous example to the integers is the 
QRW associated with the unitary operator 



(27) s = £ \i + ® ITXTI + E I* - ® U>UI- 



V 



/ 
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The evolution is now even simpler than the previous one: all spins up 
move to the right and all spins down move to the left. 
Ordering the pure states as follows 

• ••, I - 1) ® IT), I - 1) ® II), |o> ® IT), |o) ® ||), |i) ® IT), |i) ® II),- • • 

the transition matrix U is the doubly infinite matrix 












1 














1 






































1 














1 






































1 














1 






































1 














1 












Clearly the model is more complicated than it needs to be: we can 
separate it into two non-interacting evolutions corresponding to spins 
up and down. Consider for instance the set of spins pointing up which 
move in one unit of time one step to the right. The corresponding 
matrix is now the doubly infinite matrix 

1 
1 

0_J 

1 
1 
1 

The doubly infinite column vector v (z) = (..., z~ 2 , z~ l , 1, z, z 2 , . . . ) T 
clearly satisfies 

U^v{z) = zv{z). 

Since the components Vj(z) = z\ j G Z, of v(z) are orthonormal with 
respect to the Lebesgue measure, arguments entirely similar to those 
of the last paragraph of the previous example show that 

r d z 

v k (z)^vj(z)— = (£?);,*, 

JJ Z,HIZ 

an identity which can be trivially checked by a direct computation. 
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There is another approach to the problem that reduce this example 
to the previous one. Since all spins are pointing up we can ignore their 
orientation and concentrate on their locations . . . , —2, —1,0, 1,2,... as 
a way of describing the pure states. There is now a useful trick (used 
in a previous example in Section [2]) that consists in "folding" the set 
of all integers in the way 

0,-1,1,-2,2,-3,3,... 

or equivalently relabelling them as follows 

...,5,3,1,0,2,4,6,... 

With this relabelling the doubly infinite matrix U-\ becomes the semi- 
infinite matrix U + and we are back in the situation discussed in the 
previous example. 

The evolution of the states with spins pointing down is also de- 
scribed by the matrix U + of the previous example, but the required 
"folding" in this case is 

-1,0,-2,1,-3,2,-4,... 

Notice that the initial transition matrix U is the doubly infinite 
CMV matrix with null Verblunsky parameters. We will see that com- 
mon QRWs on the integers, such as the Hadamard one, essentially 
correspond to other less trivial doubly infinite CMV matrices. 

The elementary character of the previous examples is the reason 
to introduce them early on, since they show in an extremely simple 
way the typical features of the method that we will employ in more 
intrincate cases. For instance, the decoupling of a QRW on the integers 
into two QRWs on the non-negative integers will take place in the 
next examples too, a fact that simplifies considerably their analysis. 
Nevertheless, the decoupling in the following examples is not so easy 
to notice since it holds in a basis of mixed states instead of pure ones. 

7. A MORE GENERAL CASE: QRWS WITH DISTINCT COINS 

Here and in the rest of the paper "distinct coins" means not neces- 
sarily "constant coins". 

The previous example on the integers can be generalized to more 
interesting QRWs with no decoupling between up and down states, by 
including possible transitions between such states. A simple way to do 
this is to consider the following dynamics: a spin up can move to the 
right and remain up or (finally we get away from deterministic models) 
move to the left and change orientation. A spin down can either go to 
the right and change orientation or go to the left and remain down. 
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In other words, only the nearest neighbour transitions such that the 
final spin (up/down) agrees with the direction of motion (right/left) are 
allowed. This dynamics b 63XS cL resemblance to the effect of a magnetic 
interaction on quantum system with spin: the spin decides the direction 
of motion. 

Schematically, the allowed transitions are 



l<>®lt> 

l«')®li> 



\i + 1) ® |f) with probability amplitude c l n 

\i — 1) <8> |4) with probability amplitude c l 21 

\i + 1) ® |t) with probability amplitude c\ 2 

\i — 1) <g> IJ.) with probability amplitude c l 22 



where, for each i 6 Z, 



V C 21 C 22/ 



is an arbitrary unitary matrix which we will call the i coin. 

Notice that this is already more general than the Hadamard exam- 
ple usually discussed in the literature. We will analyze the general case 
of a constant quantum coin in more detail in the next section, and then 
the popular Hadamard special case in a later section. 

The transition matrix U is the unitary doubly infinite matrix 



-21 
-,-2 
•22 













21 

-22 



-11 
.-2 
12 












c 21 
L 22 



11 
12 



70- 

c ll 

c 12 















C 21 
C 22 









-11 
,1 

12 



22 



CGMV 



which has the structure of a doubly infinite CMV matrix with null odd 
Verblunsky parameters 



U Ct_4 


n 

U 


n 
u 


P-4 














P-4 








— a_4 




















a -2 








P-2 














P-2 








— a_ 2 




















a 








Po 














Po 








-a 




















«2 








P2 














P2 








-a 2 



However, U is not exactly a CMV matrix unless the diagonal elements 
of every coin Cj are positive so that they can be identified with p 2i . 

Notice that |c^| = |c 22 | and \c\ 2 \ = |c 21 | due to the unitarity of Cj. 
The coin Ci will be called trivial when c* n = c 22 = and non trivial 
otherwise. If some coin is trivial then the transition matrix U becomes 
a direct sum of two unitary matrices. In what follows we will assume 
that every coin is non trivial, otherwise the QRW splits directly into 
independent simpler ones in the basis of pure states. 

Now, a simple change of phases in the basis transforms U in a CMV 
matrix: if e*°fc is the phase of c? kk then it is easy to check that C = A^UA 
is a doubly infinite CMV matrix where A = diag(. . . , A_i, Ao, Ai, . . . ) 

is given by A 2 j+ 2 = e - * "^^ and A 2 j+i = e* cr 2A 2 j-i with A_i = A = 1. 
The corresponding Verblunsky parameters are 

«2j = C 2 1\ , «2j+l = 0, 

A2j-1 

so that p2j = |cii| are positive. 

As in the previous example on all the integers, we can do a folding 
trick to transform the doubly infinite transition matrix in a semi-infinite 
one. The appropriate ordering of the pure states to obtain a banded 
semi-infinite matrix as narrow as possible is 

|o)®|T),|-i)®U),|-i)®IT),|o)®U>, 

(29) |1>®|T),|-2)®U),|-2)®|T),|1)®U), 
which combines in a suitable way the orderings 

|o> <s> |T>, I — 1> <s> |T>, |i> <s> IT>, I — 2> <s> |T>, . . . 
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and 



1>®|I>,|0)(8)||),|-2)®||),|1)®||),... 



used in the example discussed in Section The reason for this choice 
is that, in contrast to that case, up and down states do not decouple 
now and we must interlace their orderings. 

Let us denote by U, C, A the result of performing such a reordering 
on U,C,A. Then, C = A^UA, with A diagonal unitary and U,C with 
a 2 x 2-block CMV structure. Moreover, C is exactly such a block CMV 
matrix, and its matrix Verblunsky parameters ( oij)°^_Q are 



a 



2j 





a 2j 



-Oi-2j-2 





OL 



2j+l 



o. 



where (aj)j G 2 



are the Verblunsky parameters of C. Explicitly, 



(30) 








P L o 








Po 





-a Q 












ol\ 








P L 2 






P2 








Ot-2 












ol\ 









V 



Pa 



P\ 

-a 4 



7 



The positive definite matrices pf-, p^,- are given by 



pi 









P2j 



pi 



P2j 







p2j-2 



Summarizing, the class of QRWs on the integers with arbitrary non 
trivial quantum coins can be described either by using doubly infinite 
CMV matrices or 2 x 2-block CMV matrices, in both cases the odd 
Verblunsky parameters vanish. The interest of these results is that 
they allow us to obtain for some examples of these QRWs a KMcG 
formula in terms of matrix valued Szego polynomials. 

Such a KMcG formula comes from the fact that the block CMV ma- 
trix C has an associated vector x = (1, x\, a^, . . . ) T of left orthonor- 
mal Laurent polynomials satisfying Cx(z) = zx(z). Let Aj be the 
2x2 diagonal blocks of A = diag(A , A 1; . . . ). Then, A = 1 and 
UX(z) = zX(z) for X = (1,X 1 ,X 2 ,...) T defined by X j = A jXj . 
Besides, {Xj)°^ is left orthonormal with respect to the orthogonality 

oo 



measure p of (x 



3 > j=o 



i.e., 



{Xj)jt and (xj)^L are left orthonormal 
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Laurent polynomials for the same measure with different normaliza- 
tions). In consequence, 

» OO „ 

/ z n X j {z)d^{z)X k {z)^ = E / (U n )j,iXi(z)d^(z)X k ( Z y = (U n )j, k , 

1=0 

which can be interpreted as a KMcG formula for the related QRW 
on the integers. Here {U n )j k stands for the 2x2 blocks making up 

(U n )n- n =( U ^ M U lh-(k+l)l 
\U )2j,2k I n „.n 

(31) 



2j,2fc+l - \ n ..n 

' (j+i)4,-(fc+i)T -U+VIM 

'Tjn\^ _ ( U -{j+l^M U -(j+l)T ,-(fc+l)J, 

i u J2j+l,2fe — I 



u il-{k+i)] u jlM 



u™ k being the probability amplitude to go from the pure state j to the 
pure state k in n steps. 

The measure fx and the orthonormal Laurent polynomials {X j)°?Q 
in the above KMcG formula will be called the measure and orthonormal 
Laurent polynomials associated with the related QRW on the integers. 

8. QRWS WITH A CONSTANT COIN 

The Hadamard QRW is an example of the QRWs described in the 
previous section. It corresponds to a constant coin Cj = H given by 

1 A 1 



(32) H -fiV -i 

The Hadamard QRW is an example of an unbiased QRW, i.e., a QRW 
with a constant coin such that all the allowed transitions are equiprob- 
able. 

We will consider in this section a more general class of QRWs on 
the integers: those which have an arbitrary constant (unitary) coin 



(33) a = C 



C\l Cx2 

Identifying C with the operator on the spin state space given by 
(34) C|t> = cult) + C21II), C\i) = c 12 |t> + csoli), 

the QRW evolution operator can be written S <E> C with S as in ([2 
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Let e la ~ k be the phase of c^. According to the previous discussion, 
the doubly infinite transition matrix is U — ACA^, where 

A = diag(...,A_ 1 ,A ,A 1 ,...), \ 2j = e~ ijr7 \ \ 2j ^ = e l ^\ 
and C is the doubly infinite CMV matrix with Verblunsky parameters 

a 2j -i = 0, a 2j = c 21 e~ l2j ®, = -(d + a 2 ). 

The semi-infinite form U of this transition matrix can be factorized 
as U = AC A" 1 " with A = diag(l, A 1; A 2 , . . . ) given by 

and C the 2 x 2-block CMV matrix with Verblunsky parameters otj = 
A je -^ +1 ^, where 

(o odd j, . f0 -a\ _ w 

J [A even j, V a / 

The Verblunsky parameters ctj are simultaneously unitarily diago- 
nalizable because 

Hence, we can use the infinite matrix P = diag(P, P, P, . . . ) to trans- 
form C into a new block CMV matrix C (0) = P j CP with diagonal 
Verblunsky parameters given by 




where 

(35) cij 



odd j, 
i\a\ even j. 



This means that, although the pure states do not decouple, there 
is a basis in which the QRW decouples in two independent ones corre- 
sponding to scalar CMV matrices with Verblunsky parameters 



More precisely, let and fi± be the orthonormal Laurent polyno- 

(o) 



mials and measure associated with the Verblunsky parameters ±a 
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Since U = QC (0) Q j , Q = AP, we find that UX(z) = zX(z) for 
X = (1, Xi, X 2 , . . . ) T defined by 

and Xj are left orthonormal with respect to the same matrix measure 
H as Xj, that is, 

" - p ("o + **■ 

Therefore, 



)j,ki 



(36) / z n X j (z)d f j,(z)X k ( Z y = (U r ' 
Jj 

which will yield a KMcG formula for the QRW with a constant coin on 
the integers once the measure \x and the orthonormal Laurent polyno- 
mials Xj are obtained in Section [TUl 

So far, we have split a QRW on the integers with a constant coin 
into two QRWs on the non-negative integers. Now, to obtain a KMcG 
formula we just need to compute the scalar orthonormal Laurent poly- 
nomials and the corresponding scalar orthogonality measures fi±. 

Moreover, the corresponding Verblunsky coefficients are ±a 7 -e" -t ^' +1 ' 1 ' 
with a,j given in (135]) . hence (fl2|) shows that [x± can be obtained rotat- 
ing by an angle $ the measure with Verblunsky parameters 

(37) ±z|a|, 0, ±z|a|, 0, ±z|a|, 0, . . . (|a| = |c2i|) 
Explicitly, from (1121) and (1131) we find that 

d/*(s)=dA(e-"s), » = P (^o A °_) pt > 

(38) Xj{z)=k j x j {e- i »z) 1 *i = p(f 

/ e y(o-i-o- 2 )/2 q \ „ A 

A 2i- 1 = ( o e i - ?(<T2 - CTl)/2 J ' A2i = A ^'- 1 ' 

where xf and /t-t are the orthonormal Laurent polynomials and mea- 
sures with Verblunsky parameters (1371) . Hence, the matrix Caratheodory 
function F of /x is related to the Caratheodory functions F± of fi± by 

(39) F(z)=P(e-»z), F = p[ F q p)^- 

A CMV matrix with the Verblunsky parameters (|37j) can be un- 
derstood as the transition matrix of a QRW with a constant coin on 
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the non-negative integers. Thus, a QRW with a constant coin on the 
integers always splits in two QRWs with constant coins on the non- 
negative integers. Such a splitting is not initially obvious since it takes 
place in a basis of mixed states. 

It is worth to remark that the block CMV matrix related to a 
constant coin C on the integers depends only on the modulus of the 
entries of C (one of them determines the others by unitarity). This 
means that QRWs on the integers with a constant coin do not depend 
essentially on the phases of its entries. For an unbiased QRW all the 
entries of the coin C have equal modulus, so \a\ = \c 2 i\ = l/y/2 due to 
the unitarity of C. The Hadamard QRW is therefore just a canonical 
example of unbiased QRW on the integers. 

Bearing in mind that we are reducing the problem on the integers to 
two problems on the non-negative integers, in the next section we will 
discuss completely the QRWs with a constant coin on the non-negative 
integers. 



9. QRWS WITH A CONSTANT COIN ON THE NON-NEGATIVE 

INTEGERS 

QRWs with a constant coin on the non-negative integers are a pre- 
liminary step to complete the discussion of QRWs on the integers. Be- 
sides, these QRWs have their own interest because they present special 
features that do not appear for QRWs on the integers. 

Let us suppose that we order the pure states on the non-negative 
integers as in (1261) . A unitary matrix like 



C 21 













°22 





c 12 












C 21 








l 1 






C 22 








C 12 












c 2 

c 21 

r 2 

c 22 












v ■•• ■•• ■•• ••• ■■■/ 

can be understood as the transition matrix for a QRW on the non- 
negative integers with arbitrary (unitary) coins Cj as in f[2"gj) for i = 
0,1,2,.... 

A constant coin C corresponds to the operator S + ® C where S + is 
given in (1251) and C is identified with the operator (1M1) . 

Exactly as in to Section [3, even in the case of (non trivial) distinct 
coins, a simple change of phases in the basis reduces U to a scalar CMV 
matrix with null odd Verblunsky coefficients: C = A^UA is CMV for 



28 



CGMV 



A = diag(A , Ai, . . . ) given by \ 2j+ 2 = e l0 iA 2 j, A 2 j+i = e^X^j-i with 
A-i = Ao = 1 and e io fc the phase of c? kk . Its Verblunsky parameters are 
given again by 



A 2j-1 

so that p 2 j = |cii| are positive. 



_ -j *2j _ „ 

Q 2j — C21T > «2j+l — u, 



If /x and Xj are the measure and orthonormal Laurent polynomials 
associated with C, then UX(z) = zX(z) with X = (1, X 1; X 2 , . . . ) T 
and Xj = XjXj. Since /x makes the polynomials Xj orthonormal too, 



which is a KMcG formula for a QRW with arbitrary distinct coins on 
the non-negative integers. 

The above formula simply says that the transition matrix U has 
been identified as a matrix representation of the unitary multiplication 
operator on L 2 (Y) given by 

The measure /1 and the orthonormal Laurent polynomials Xj in 
(140]) will be called the measure and orthonormal Laurent polynomials 
associated with the related QRW on the non-negative integers. 

Following the arguments and notation of Section [8] we can see that 
a constant coin C leads to a rotation of the case associated with a CMV 
matrix with Verblunsky parameters 



In contrast to (j57|) . the even Verblunsky parameters in (1411) can have 
an arbitrary phase which comes from the phases of the coefficients of 
the coin. Hence, we should expect a stronger influence of the phases 
of the coin for constant coins on the non-negative integers compared 
to the same situation on the integers. This will bring in new possibili- 
ties when discussing the semi-infinite version of the Hadamard QRW, 
corresponding to \a\ = l/v2. 

To be more precise, if A = diag(Ao, Ai, . . . ), A 2 j = e _y<T1 , A 2 j-i = 
e u<72 , then C = AWA is a CMV matrix whose measure \x is obtained by 
rotating by an angle d the measure /2 with Verblunsky parameters (1411) . 
If Xj and Xj are the orthonormal Laurent polynomials associated with /x 
and /x respectively, we know from (fT3|) that X2j~i(z) = e~ lJ X2j-i(e~ z) 
and X2j(z) = e t ^X2j(e~ t ^z). So, according to the previous discussion, 
the Laurent polynomials 



(42) Xj(z) = XjXjiz) = Xj Xjie-^z), X 2j -i = A 2j = e^'^ 2 , 



(40) 




(41) 



a, 0, a, 0, a, 0, . . . 



a = c 2 ie 
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are orthonormal with respect to /i, satisfy UX(z) = zX(z), X = 
(1, Xi, X2, . . . ) T , and provide the KMcG formula (|40p . which can be 
rewritten as 

(43) (U n ) hk = e in ^ [ z n x 3 (z)J^z)dfi(z). 



To get from ( 1401) a KMcG formula for the QRWs with a constant 
coin on the non- negative integers we simply need Xj and ft, which are 
calculated in an appendix at the end of the paper. The main results 
from the appendix are summarized now: 

An expression for the orthonormal Laurent polynomials is given by 

x 2j -i(z) = Uj{y) - p- l {z + a)Uj- X {y), _ 1 

x 2j (z) = U j (y)-p- 1 ( Z - 1 + a)U j _ 1 (y), V ~2p [Z Z h 



with p = a/1 — |a| 2 and Uj the second kind Chebyshev polynomials 
given by (j4j. In particular, x 2 j(z) = x 2 j_i(l/z). 
On the other hand, the measure is given by 

d9 

dp(e id ) = w(9) — + M5(6 - p)d6, 

Z7T 



v^in 2 9 — sin 2 77 . . . r , 

W ^) = ^T— Q (9 G 77,7T- T? U 77 -7T, -77 , 

'suit/ - snip I 



M 



I Rea| a/ sin 2 77 — sin 2 (3 



a/1 - |Ima| 2 I cos/9| 

where the angles 77 G [0, tt/2) and /3 G (— 7r, tt] are defined by means of 
sin 17 = I a I , sin/3 = — Ima, sign(cos/3) = sign(Rea). 
The corresponding Caratheodory function is F(z) = F(e~ l ®z) where 



p sj(z - z- 1 ) 2 + 4|a| 2 + 2Rea _ z - z" 1 - 2ilma 

r [Z) = 



z - z- 1 + 2ilma ^/( z - z^ 1 ) 2 + 4|a| 2 - 2Rea 

Now that we have all the ingredients that enter in the integral in 
l3|) . notice that the only parameter that appears in it is the value 
of the complex number a. The parameter d as well as the individual 
values of <J\,<J 2 have an effect on {U n )j t k but this appears only as a 
factor in front of the integral on the right hand side. 

The measure d/i(z) = d\x{e~ %{> z) appearing in (J4D]1 has in general 
a continuous weight plus a Dirac delta. The weight is supported on 
two symmetric arcs of angular amplitude 2r] centered at ±ie l1? . The 
mass point is located outside of the support of the weight and it is 
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absent only when a is imaginary, which holds exactly when the coin is 
symmetric (c 12 = c 2 i). 

While the location of the weight only depends on \a\ = |c2i|, its form 
and the location and mass of the Dirac delta also depend on the phase 
of a = C2ie . Therefore, the phases of the coin C have a remarkable 
influence in the semi-infinite QRWs with a constant coin, in marked 
contrast to the case of the integers. 



10. QRWS WITH A CONSTANT COIN ON THE INTEGERS 

The results of the previous section permit us to complete the analy- 
sis of the QRWs with a constant coin on the integers. Let us remember 
the notation: C = (0^)^-^=1,2 is the constant coin, e lak the phase of c^, 
d = + o"2)/2 the angle of rotation, a = ~C2\e l ® and p = \J\ — \a\ 2 . 

The matrix ingredients Xj, n of the KMcG formula fl36|) are given 
by (|38l) in terms of the scalar ones x^, fi± with Verblunsky parame- 
ters ±i|a|, 0, ±z|a|, 0, . . . , a specialization for a imaginary of the case 
a, 0, a, 0, . . . analyzed in the appendix. This is precisely the case 
where the mass point disappears, so we can anticipate that the matrix 
measure for a constant coin on the integers is always given exclusively 
by a continuous matrix weight. 

Combining the results of Section [8] and the appendix we get the 
following expressions for the scalar objects 



- y/(z - z- 1 ) 2 + 4|a| 2 z-z~ l ^2i\a\ 

b±\z) - - 



z-z- 1 ±2i\a\ y(^-^- 1 ) 2 + 4|a| 2 ' 



,~ / ids ,„sd9 / sin 9 =p sin n 

2tt y sinffi smr] 

9 E [t], 7r — 77] U [7] — 7T, —7]}, smr] = \a\, T] E [0, 7r/2), 



x 



4 (z) = UM- p-\z ± i\a\)U^{y), = j_ 
x%(z) = x 2j -i(l/z), y 2p 

so that the matrix objects in (1381) and ( 1391) are given by 

F(z) = F(e-^z), 

F(z) ^ fz ~ Z ~ l 271 



y/(z - z~ 1 ) 2 + 4|a| 



i\2 \ A\ n \2 \ —2a z — z 
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lix 





d£i(e ie 


sin 9 




±ia 


sin 9 J 



' 9 E [q- 71,-7]}, 



Xj{z) = A jXj (e-^z), A 2i _! = i^ 3{ai ^ )l2 eijia2 °_ ai)/2 \ = At., 
i 2 i-i(z) = lUj(y) - ^ ^ f^jUj.^y), x 2j (z) = x 2j -i(l/zy . 

The expression for the polynomials Xj can be compared to the one 
at the end of Section [2j 



11. The Hadamard QRW versus other unbiased QRWs 

The Hadamard QRW is the unbiased QRW on the integers with 
constant coin ( |32l) . Applying the results of the previous section to this 
case gives 

1 fl + z> V2z 



VTTz 1 \V2z l + z' 



2 / ' 



7T 7T1 



iQ\_ 1 /"VI + cos 20 ±1 \ d0 . f J^eL-4'4J' 



= X 2j _y{l/Z), 

where, here and below, 



If 



z — z 1 



V2i 



and the square root y/1 + is the analytic branch with value 1 at the 
origin. Notice that the measure /x is symmetric with respect to the 
real line and the matrix coefficients of the Laurent polynomials Xj are 
real. 

We can compare this case with a semi-infinite version with the same 
constant coin (!32|) . Taking into account Section [9j the coin (1321) yields 
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on the non-negative integers 



F(z) 



1 - V2 



z + z 



2 



, , m v / l + cos2^ ±1 d9 
V cos 26* 27r 




= i J {Z7j(j/*) + t(V2z - 1)^-1(2/*)}, 
X 2i (z) = (-l)^_!(l/z). 

We can observe the similarity between the Hadamard QRW and its 
version on the non-negative integers. Again, the measure \i is given by 
a symmetric weight and the Laurent polynomials Xj have real coeffi- 
cients. 

However, this similarity goes away for other unbiased QRWs. Con- 
sider for instance the constant equiprobable coin 

(44) c= ^G : i 

Section [10] gives for the corresponding QRW on the integers 
1 + z 2 iy/2z 



F z) 



VTT^ \-iV2z 1 



z 




, i9) _ 1 / VI + cos 29 ±i \ d6_ 
X % -(z)=X^_ 1 (l/^). 

As in the Hadamard QRW, the measure /x is given by a weight which 
is symmetric with respect to the real line. 

From Section [9j the same constant coin (j44p yields for the non- 
negative integers 

F(z) = - 



l + z 2 ^/TTz* + iV2z' 



cos 20 d# 1 
+ cos 29 2tt + ~j 

X2j-x(z) = t'iUjiy*) + i{V2z + ip^iy*)} 



Me*) = d r ^- ^ + T2 S(0 ~ - 2 )d0, e h| , | ] U ft, f ], 



x 2j ( z ) = i-iyx^ii/z). 
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The differences between the infinite and the semi-infinite versions 
of this unbiased QRW are evident. On the non- negative integers the 
weight is symmetric with respect to the real line, but there is also a 
mass point which destroys such a symmetry for the full measure. 

The similarities and differences observed above regarding the ingre- 
dients of the KMcG formulas for different QRWs reveal how close or 
far their probabilistic behaviours are from each other. To make this 
more evident let us obtain the probability amplitudes u^ k for a n-step 
transition between certain states j, k in the different examples. 

For the examples on the integers, according to (1311) . 

(TT n \ — ( u ot,ot m ot ,-U 
\ u -ii,0T "-U-U 

Due to the translation invariance of any QRW with a constant coin 
on the integers, (C/ n ) ,o gives the probability amplitudes 11%+^, u k , k ^, 
< T , (fc _ m and < Ufc+ i )T for any k. 

On the other hand, the KMcG formula states that (U n )o,o = 
J T z n dfj,(z) = /-i n are the moments of the related measure, which, fol- 
lowing f)24p . are provided by the McLaurin series of the Caratheodory 
function F. 

For the Hadamard QRW on the integers the expression of F gives, 
apart from the trivial moment fi = 1, 

- Cm i /0 A 

/*4m ~~ M4m+2 ~~ ^ ' A*4m+1 ~~ \\ Q J ' / X 4m+3 ~~ U, 

where 

(45) c = l, ^ = (-l) n Tlfl-iV n>l, 



k=l 



are the coefficients of the series 1/yT+z = ^2^L c n z n . Hence, 



n=4m,4m+2, 
n=2m+l, 



I ~rs n=4m+l, 
7/ n — n n — I v2 



n=2m,4m+3, 



The non-diagonal elements of 



n oi,ot w oj.-u 



permit us to complete the description of the n-step transition ampli- 
tudes between the spin states at the same site. From (}36l) and the 
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expression of the polynomials X k for the Hadamard QRW we get 



i o) Mn + v ^(o -1 



thus, the values of the Hadamard moments give 



%T,(fc+l)T ~~ U kl,{k-l)l 




n=2m,4m+l, 
n=4m+3, 



except for u£ Tj(fc+1)T = -Uki,(*-i)i = Vv^- 

In the Hadamard example, all the n-step transitions between pure 
spin states for the same site are forbidden for odd n (except for n — 1), 
while for even n we find probabilities |c m | 2 /4 if n = 4m, 4m + 2. The 
same result holds for any other unbiased QRW on the integers due to 
the similarity between the corresponding Caratheodory functions and 
orthonormal Laurent polynomials. 

Consider now the Hadamard coin on the non-negative integers. 
From the McLaurin series of the associated Caratheodory function and 
(113]) we find that, apart from fi — 1, the related moments are 



_ c m ^ c m ^ 



The KMcG formula (1401) and the expression X x (z) = 1-V2Z- 1 for 
the first orthonormal Laurent polynomial yield 



(^™)o,0 = (J>n, (^ n )l,l = 3jU n - V2fl n +1 - V^/in-l, 
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Therefore, 



ot,ot 



u 



ot.oj. 



Cm 

— n 


=4m,4m+2, 


Cm 




=4m+l, 
=4m+3, 


Cm. 

2 


n=4m, 





n=4m+l, 


Cm 

2 

c m +i 


n=4m+2, 
n=4m+3, 



01,01 



It 



oi,ot 



Cm 

— 

Cm 

72 

Cm+C 7n 






C ; 1 1 

2 


n=4m, 





n=4m+l, 


Cm 

2 


n=4m+2, 


C??i 


n=4m+3, 



except for u\ 



and «oj,,ot 



1/V2. 



The above results show that, concerning the n-step transitions be- 
tween pure spin states at site 0, the Hadamard coin on the non-negative 
integers has the same probability amplitudes as the Hadamard coin on 
the integers for even n, while for odd n some of the transitions remain 
forbidden. 

On the other hand, the Caratheodory function for the unbiased 
QRW on the non-negative integers with coin (144)) yields, apart from 
Ho = 1, the moments 



d 



A*-4r 



m 

2~' 



A*4m+1 — 7y^' / i 4m+2 



d 



in 
2~ : 



/Mm+3 



where d n are the coefficients of y/1 + z/(l — z) = J2n°=o dnZ n , so that 
V / TT^/(1 + z) = Y.n=Mz 2n - z 2n+l ). Explicitly, 



(46) 



d n = c + ci H h c n , VT+z = 

f 1 \n-l n 



n=0 



c — 1) Ci — — ■, 



k=2 



n > 2. 



Introducing -Xi(z) = — (i + V^-z x ) in the KMcG formula fT40l) gives 



(^ n )o,o = A*n, (^ n )i,i = 3 /i„ + iV2fi n+1 - zy^^n-i, 

(^™)l,0 = ~(^n + V^/in-l), (f /n )o,l = ¥n ~ V^/Vf-l- 
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Hence, 



( d„ 



U 



ot,ot 



it 



01,01 



V2 



—i 



v/2 

2- 



n=4m, 
n=4m+l, 
n=4m+2, 
n=4m+3, 

a n=4m, 



/ 3d„ 



It 



OJ.,01 



2 

,- 3-2rf m 
' V2 

3d m 

2 

; 3— dn 



—l- 



V2 



-1 



' 2 



n=4m+l, 
n=4m+2, 
n=4m+3, 



oj.,ot 



V2 
;2-d„ 



n=4m, 
n=4m+l, 
n=4m+2, 
n=4m+3, 

n=4m, 
n=4m+l, 
n=4m+2, 
n=4m+3, 



except for wJj.,oj = and Mqj,,ot = — Vv2- 

The probability amplitudes of this example are quite different from 
those of the Hadamard coin on the non-negative integers. They also 
show quite a different behaviour when compared to any unbiased QRW 
on the integers, including the case of the same coin 044]) . In particular, 
if n > 2, this example has no forbidden n-step transitions between the 
pure spin states at site 0. This is due to the inequality (I2 < d n < d\ 
for n > 3, which is a consequence of (146]) and the fact that \c n \ is 
decreasing and c n has alternating signs for n > 1. 



12. Asymptotics of QRWs 

In Sections [7] and [9] we have seen that the transition matrix of any 
QRW (on Z > or Z) with non trivial distinct coins has an asso- 
ciated (scalar or 2 x 2-matrix valued) measure. As in the previous 
section, the corresponding Caratheodory function allows us to com- 
pute the moments of the measure which, with the aid of the related 
orthonormal Laurent polynomials, provide through the KMcG formula 
the amplitudes of the n-step transitions for any value of n. This opens 
the possibility of studying the asymptotic behaviour of such amplitudes 
when n goes to infinity. 

Several different authors have obtained specific asymptotic results, 
mainly in the case of the Hadamard walk, by using different methods. 
For a very good account, see [Tj5]. We have not attempted any compar- 
ison between our rather general results and the many detailed results 
in the literature. The results in this section are given to indicate how 
our method could be used for similar purposes. 

For instance, in the case of the Hadamard coin on the non-negative 
integers, as well as for any unbiased QRW on the integers, the n th 
moment has zero limit when n goes to infinity because the coefficients 
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c n in (H5l) satisfy lim^oo c n = 0. Therefore, lim^oo ^ = in these 
cases. 

On the other hand, the unbiased QRW on the non-negative inte- 
gers with coin (jUJ) gives moments /i n satisfying linv^^ \fj, n \ = l/y/2 
because, from (HoT) . 



lim d n = } c k = a/1 + z\ = \/2. 

k=0 

Hence, the probabilities for returning to the spin up 0-state, if the sys- 
tem was originally in such a state, converge to a non-zero limit, i.e., 
lim^oo Iuq^otI 2 = 1/2- Indeed, although the probability amplitudes 
m ot,ot = A*n do not converge in this example, the quantities i~ n \i n actu- 
ally converge to 1 / a/2 when n goes to infinity. 

The moments are given by only a few of the coefficients of the pow- 
ers of the transition matrix U: the coefficient (0, 0) of U n is the scalar 
moment \x n in the case of the non-negative integers, while the coeffi- 
cients (0,0), (0,-1), (—1,0) and (—1,-1) of U n provide the matrix 
moment n n for a QRW on the integers. However, as we saw in the pre- 
vious section, the rest of the transition amplitudes can be calculated 
in terms of the moments using the KMcG formulas. In fact, according 
to Sections [7] and [HI this should be possible, not only for QRWs with 
a constant coin, but for any QRW with non trivial distinct coins. As 
a consequence, for all these kinds of QRWs, the asymptotic behaviour 
of the moments controls the asymptotic behaviour of the powers of 
the transition matrix. This is the idea behind the result given in the 
following proposition. 

To get a better understanding of what is coming, we start with 
some remarks. Looking at our examples, we have to deal with the 
situation where the n th moment multiplied by some phase e~ ldn is con- 
vergent. As we will see, this condition is connected with the situa- 
tion where e~ ie "(U n )ij is convergent for any Since the sequence 
U n is uniformly bounded because U n is unitary for any n, the exis- 
tence of lim^oo e~ l9n {U n )ij for any i,j is equivalent to saying that 
lim^oo il)e~ t9n U n rf exists for any row vectors xjj,r) G L 2 (I x S), which 
defines the familiar weak convergence of operators, see [T71 Chapter 
III]. 

The weak limit U°° = (C/°°), U°° = lim^oo e- i6n (U n )ij, when it 
exists, defines an operator on L 2 (I x S) which provides the asymptotic 
behaviour for n going to infinity of the n-step transition amplitude 
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between any two states i/j,r) because 

lim UU n rf - e ie "^°y) = 0. 

n— >oo 

Although U°° does not inherit in general the unitarity of U n , it has a 
norm not greater than one because U n does so for any n. 

Proposition 12.1. Let U be the transition matrix of a QRW on the 

integers or the non-negative integers with arbitrary non trivial distinct 
coins. Concerning the asymptotic behaviour as n goes to infinity of U n 
and the n th moment of the related orthogonality measure, we have the 
following results: 

(1) U n converges weakly to zero if and only if the n th moment con- 
verges to zero. 

(2) For any sequence e ldn of phases, e~ t9n U n has a non null weak 
limit if and only if the n th moment multiplied by e~ l0n converges 
to a non null limit and lim^oo e^ 6n+1 ~ 9n ^ exists. 

Proof. Consider a QRW on the integers. Let U be the 2 x 2-block five- 
diagonal matrix obtained by performing the folding ( 1291) on U. Then, 
U n converges weakly to zero if and only if U n does so. 

The KMcG formula obtained in Section [7] states that the 2x2- 
blocks of U n are given by 

{U n ) WJ = [ z n X l (z)d f i(z)X ] (z)\ 
Jt 

where X{ are the corresponding orthonormal Laurent polynomials. If 
Xi(z) = Ylt—p Ai±z k with 2x2 matrix coefficients A^, then 

(U n )i,j = ^ ^ A it ktJ- n+ k-iA\ d . 

k=pi l=pj 

This equality implies that lim n ^, OQ (U n )ij = for any i,j exactly when 
hm n ^ oc n n = lim n _ >00 (C/ n ) ,o = 0, which proves (1). 

Concerning (2), notice that e~ l9n U n and e~ %dn U n have a non null 
weak limit simultaneously. Suppose that both lim^oo e~ ldn ii n and 
lim^oo e l ( 9n+1 ~ 9n ^ exist. The above equality gives 

k=Pi l=Pj 

If lim n ^oo e i ^ n+1 ~ en " > = z , then hm n ^ 0O e l ^ n+k ~ 6n ^ = z*. Thus, the 
previous identity shows that \\xa n ^ 00 e~ ldn {U n ) i j exists for any i,j. 
Moreover, if lim^oo e~ i9n n n ^ 0, then lim n ^oo e~ idn (U n )ij ^ at 
least for % = j — 0. 
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Conversely, assume that lim^oo e~ i9n (U n )ij exists for any i,j and 
is non null for some Then, lim^oo e~ ldn n n = lim n ^oo e~ t9n (U n )ofi 
exists and must be non null because otherwise lim^oo e~ l9n (U n ) it j = 
for any i, j due to (1). Denote by U°° the (non null) weak limit of 
e -i&njjn T a ki n g W eak limits in Ue~ idn U n = e i (^+i-^) e _i0n + 1 C/ n+1 we 
obtain UU°° = \im n ^ 00 e t( - en+1 ^ e '^U 00 , thus e^ en+1 ~ dn] must converge 
because U°° is not null. 

The proof for a QRW on the non-negative integers is completely 
analogous and even simpler because we do not need the folding. □ 

As a consequence of the previous results, the transition matrix of 
any unbiased QRW on the integers converges weakly to zero, which 
means that the amplitude of the n-step transition between any two 
(finite or infinite) superposition of pure states converges to zero as n 
goes to infinity. This is also true for the Hadamard coin on the non- 
negative integers. On the contrary, the transition matrix U for the 
unbiased QRW on the non-negative integers with coin 044|) should be 
such that i~ n U n converges weakly to some non vanishing weak limit 
U°°. Indeed, we can compute such a weak limit with the aid of the 
following result. 

Proposition 12.2. Let U be the transition matrix of a QRW on the 

non-negative integers with arbitrary non trivial distinct coins. If, for 
a sequence e tdn of phases, e n U n has a non null weak limit U°° , the 
related orthogonality measure /i has a mass point z to which e l ( en + 1 ~ dn ' > 
converges and 

U°° = fi oo X(z )X(z )\ /loo = lim e-^> n , 

n— *oo 

with X = (1, Xi, X2 . . . ) T the associated column vector of orthonormal 
Laurent polynomials and \i n the moments off/,. Furthermore, 

lim P - ie "z n - 

11m e z Q — , 

rwoc MW) 

so that lim^oo ZQ n fi n = /i({z }) and z§ n XJ n converges weakly to 

v({z })X(z )X(z o y = ^— X(z )X(z )\ 

which is the orthogonal projection onto the eigenspace of U associated 
with the eigenvalue zq. 

Proof. Suppose that e~ ldn U n converges weakly to a non null limit U°° . 
From Proposition 112.11 we know that lim^oo e~ l9n fi n = fi^ 7^ and 
lim^oo e l ( 9n+1 ~ 9n ^ = zq G T. The arguments at the end of the proof of 
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such proposition yield the identity UU°° = zqU 00 . Thus, the non null 
columns of U°° must be eigenvectors of U with eigenvalue Zq. 

Let C be the CMV matrix related to U, and let x be the corre- 
sponding column vector of orthonormal Laurent polynomials. Bearing 
in mind that U = ACA) and X = Ax with A unitary diagonal, the 
comments in Section [3] show that zq must be a mass point of \i with 
a mass given by /x({z }) = 1 / 1 1 ^ ( -^o ) i 1 2 — Vll^^o)!! 2 - Moreover, the 
eigenvectors of U with eigenvalue zq must be spanned by X(zq), so the 
columns of U°° should be proportional to X(z ), i.e., 

U°° = X(z )Y, Y = (Y , Yi, . . . ) e L 2 (Z > 0). 

Notice that Y = U^ = lim™ e- ie "(U n ) , = AW 

On the other hand, e~ ien {U T ) n converges weakly to {U°°) T . The 
unitarity of U implies that z must be an eigenvalue of U T too, and 
the corresponding eigenvectors must be spanned by X(zq). Therefore, 
similar arguments to the previous ones show that 

(U°°f = Xjz^jZ, Z = (Z , Zi, . . . ) e L 2 (Z > 0), 

with Z = 

Therefore, the matrix X(zo)Z must be equal to Y T X (zq) t '. Identi- 
fying the first column of both matrices gives Y T = ZqX(zq) = fiooX^zo), 
hence U°° = /i oo X(z )X(z f ■ 

From U°° = fi OQ X( y z )X(z o y we obtain 

Voo\\X(z )\\ 4 = X(z^U°°X(z Q ) = hm X{z )^e- i9 -U n X{z ) = 

n— >oo 

= hm e-^z^\\X(z )\\ 2 , 

which proves that lim n _oo e~ l6n ZQ = ^^(zq) || 2 = /ioo/M{ 2 o})- The 
rest of the identities follow easily from this equality. □ 

The second part of Proposition 112.21 asserts that, when e~ ldn U n has 
a non null weak limit, we can suppose without loss of generality that 
e i9 n j s z n^ with. Zq the mass point of \i to which e^ 6n+1 ~ 6n ^ must converge 
and, then, /i^ becomes fi({z }). 

For instance, in the case of the unbiased QRW on the non-negative 
integers with coin (jSj), lim^oo i~ n jj n = l/v2 is the mass of the only 
mass point zq = i of the corresponding measure. According to Propo- 
sition [TO the weak limit U°° of i~ n U n is 

U°° = \ X{i)X(i)\ X = (1, X u X 2 , . . . ) T , 
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where are the associated orthonormal Laurent polynomials given in 
Section [TT] so 

X 2j _ 1 (i)=V{ U j {y/2)- (V2+l)E^_i(V2)}, 
X 2j (i) = (-iyx 2j _ 1 (i)=X 2j . 1 (i). 

Recurrence (jl]) for the Chebyshev polynomials Uj implies that Uj(\/2) = 
2 v / 2f/ i -i(v / 2) - t/ j _ 2 (v / 2), hence 

X 2j (t) = iP{(V2-l)Uj^(V2) - U 3 _ 2 {V2)} = 

= i j (V2-l){U 3 ^(V2) - (V2 + 1)^-2(^2)} = 

= z(v / 2-l)X 2j _ 2 (z). 

Therefore, X 2j -_i(i) = X 2j (z) = (i(>/2 - and 

rroo TT°° 7T°° 7T°° _ l^r) 

U 2j-l,2ifc-l — U 2j-l,2ifc _ u 2j-l,2k _ U 2j,2fc _ ^ l v z ^ 

This allows us to compute the asymptotic transition amplitude between 
any two states. In particular, denoting lim n _ >00 (a n /6 n ) = 1 by a n ~ n b n , 

;n+j-k jn+j-k 
jn+j-k— 1 jn+j— k+l 

m jt,^ ~n (v2 - ly , u nM ~ n ^_ (V2 - iy 

These asymptotic results give further indications of the different 
probabilistic behaviour of an unbiased QRW when considered on the 
non-negative integers or in all the integers. 

13. Recurrence properties of QRWs 

For a classical random walk there is an important notion that goes 
back at least to G. Polya, see [26]. We say that a state is recurrent if, 
having started there at the initial time, one returns to it with proba- 
bility one. Otherwise we say that the state is transient. For a so-called 
irreducible chain either all states are recurrent or they are all tran- 
sient. A recurrent state is called positive recurrent if the expected 
value for the time of (first) return to it is finite. When dealing with a 
birth-and-death process on the non-negative integers the correspond- 
ing orthogonality measure dm(x) has support in the interval [—1, 1] 
and plays an important role in studying these notions. 

For instance, the process is recurrent exactly when 

(47) r 1 ^i = oo. 

J -i I — x 
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Notice that this integral is the sum of all moments of the measure 
dm(x) and that the n th moment is the probability of going from the 
state to itself in n steps. 

The process returns to the origin in a finite expected time when the 
measure has a mass at x = 1. The existence of 

lim {P n ) id 

n— kso 

is equivalent to dm(x) having no mass at x = —1. If this is the case 
this limit is positive exactly when dm(x) has some mass at x = 1. 

In the classical case one gets a lot of milage out of the generating 
function S(z) of the moments of dm(x) given by 



S(z) 



dm(x) 



xz 



In particular the generating function G(z) of the sequence g n giving 
the probability of a first return to the origin in n steps 



G(z) = J2z n 9n 



n=0 



is related to S(z) by 



G{z) = 1 



S(z) 

Therefore we have that G(l) = 1 (indicating that one returns to 
state with probability one) exactly when S{1) is infinite as noticed 
above. This relation allows us to compute the expected time to return 
to state 0. This expected value is given by G'(l). 

We are confident that these ideas and results should have a nat- 
ural translation to the quantum case using the tools provided by the 
orthonormal Laurent polynomials on the unit circle. Indeed, an ap- 
proach to the asymptotics of the powers of the transition matrix using 
such a machinery has been presented in the previous section. 

We will see now that the notion of quantum recurrence can be de- 
scribed nicely in terms of the Caratheodory function introduced earlier, 
which will play in the quantum similar role to the generating 

function S(z) of the moments for a classical random walk. Neverthe- 
less, the condition for the characterization of quantum recurrence will 
be somewhat different from (|47p . 

The study of the recurrence properties in the quantum case raises 
special issues because a quantum measurement destroys the initial evo- 
lution since the system collapses into a pure state when a measurement 
is performed. Therefore, the notion of a return to a given state for the 



MATRIX VALUED POLYNOMIALS AND QUANTUM RANDOM WALKS 43 



first time in a certain number of steps has to be interpreted with care 
in the quantum case. Such an analysis in terms of a specific measure- 
ment scheme, involving an ensemble of identically prepared QRWs, has 
been proposed recently, see [30] . 

The bottom line of this analysis is that, in the quantum case, the 
recurrence of a state is characterized by the divergence of the series of 
probabilities to return to such a state in n steps. After the modifica- 
tions which are necessary to make sense of the notion of recurrence in 
the quantum case, this result is completely analogous to the classical 
one. Its importance lies on the fact that the notion of a return to 
a state in a certain number of steps is completely meaningful in the 
quantum case and its probability can be computed using the transition 
matrix of the QRW. 

More precisely, following the interpretation of the quantum recur- 
rence given in [30], a state i/j of a QRW with transition matrix U is 
recurrent exactly when 

oo 

(48) $>n(V0 = °°> p n (V) = IV^VI 2 , 

n=l 

where p n {ip) stands for the probability to return to the state ip in n 
steps. 

Consider a QRW on the non-negative integers with non trivial dis- 
tinct coins. The recurrence of the state numbered as 0, i.e., the spin 
up at site 0, is characterized by the divergence of Y^=i l(^")o,o| 2 = 
Y^=i \^n\ 2 - From the McLaurin series (1T51) of the related Caratheodory 
function F(z) we obtain 

/ \F(e^ = l + 2$> n | 2 , 

/n z — ' 

J{i n=l 

where F(e td ) = lim^x F(re l9 ), which exists for Lebesgue almost every 
9 G [0, 27r), see [271 Chapter 17]. Therefore, |0)®|T) is recurrent exactly 
when the radial limit of F(z) does not lie on L|a (T). 

The generalization of this result to an arbitrary state, as well as to 
QRWs on the integers, is the purpose of the proposition below. In what 
follows we write F(z) = F(z, dfi) when we need to make explicit the 
measure corresponding to a Caratheodory function, and similarly for 
matrix valued Caratheodory functions. Besides, we will assume that 
any scalar or matrix valued Caratheodory function is radially extended 
Lebesgue almost everywhere on the unit circle (see [5] for the matrix 
case) . 
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Proposition 13.1. Consider a QRW with non trivial distinct coins on 
1= Z or Z > 0. 

(1) If I = Z > ; let fi be the related orthogonality measure and let 
us number the states |0) ® |t),|0) <g> | ® |T),|1) t8> | !),••• as 
|0),|1),|2),|3),|4),.... ^ state |*) = EfcLo^fcl^) iS transient 
if and only if 

oo 

|/| 2 a» g L 2 m(T), f = y>***, 

2tt *■ ' 

where X k are the corresponding orthonormal Laurent polynomi- 
als. 

(2) IfT— Z ; /et /x fee £/ie related orthogonality matrix measure and 
let us number the states |0)®|T), |-1)<8>||), |-l)(g)|T), |0)®||), ...as 
|0),|1),|2),|3),|4),.... A state |#) = J2T=o^k\k) ls transient 
if and only if 

oo 

F{z,fdn?) G L|e (T), / = 

2tt ^— ' 

fe = 

where X k are the corresponding matrix orthonormal Laurent 
polynomials. 

Proof. We will only prove (2) since the proof of (1) is similar and sim- 
pler. Consider a QRW on the integers with non trivial distinct coins. 
Let U be the semi-infinite transition matrix obtained performing the 
folding (12T)|) on the doubly infinite transition matrix U. The recurrence 
of the state is equivalent to the divergence of E^L o f>nC0) , where 
p n {^) = \ijU n ip^\ 2 and ip = (-00,^!, . . .) is the wave function of 
corresponding to the same folding. 

Denoting if) k = ("02fc> "02fc+i)> the KMcG formula permits us to write 




/. uo 

= / z n f(z)dfi(z)f(z)\ f = J2^kX k . 

The above equality identifies ipU n ip^ as the n th moment of the scalar 
measure fd/j,f\ Therefore, the same arguments given before the 
proposition show that the divergence of Y^=o \4 , U n ip Jf \ 2 is equivalent 
to F(z,fdtif) i L\{T). □ 

2tt 

Notice that this proposition associates with any state a scalar func- 
tion / G L 2 (T) or a 2-dimensional vector valued function / G L 2 (T), 
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depending on whether the QRW is on Z > or Z, so that F(z, \f\ 2 dfi) 
and F(z, fdfif^) are both scalar Caratheodory functions. We will refer 
to / and / as the functions associated with the corresponding state. 

The condition (1471) for the classical recurrence at the origin depends 
only on the behaviour of generating function S(z) of the moments at 
z — 1. In contrast, the characterization of the quantum recurrence in 
terms of Caratheodory functions has to do with their global behaviour 
on the whole unit circle. 

Proposition 113.11 has the drawback that it is not given in terms 
of the Caratheodory function of the original measure associated with 
the QRW, which is the one that we directly know, but in terms of 
the Caratheodory function of a modification of such a measure. The 
following result shows that this problem can be overcome, at least 
when analyzing the recurrence of a local state, i.e., a state which is a 
superposition of a finite number of pure states. 

Notice that the function associated with any local state is a scalar 
Laurent polynomial / in the case Z > 0, or a 2-dimensional vector 
valued Laurent polynomial / in the case Z. 

Proposition 13.2. For any positive definite d x d matrix valued mea- 
sure [i and any d-dimensional vector valued Laurent polynomial f 

F(zjdfif) G Lie & f(z)F(z,d f i)f( Z y G L 2 ae. 

2tt 2ir 

Proof. If / = El= p a k z k , a k G C d , then 

F(z,fd»f) = djF^al F k (z) = I t±^t k dfj,(t). 

i,fc=P Jt Z 

Writing 

— t k = (t + z)t k ~ l + z l —t k - x 
t — z t — z 

we find that 

Fk(z) = v k + fJt k _ ± z + zF k -i(z). 

Bearing in mind that Fq(z) = F(z,d/J,), the iteration of the above 
equality yields F k (z) = L k (z) + z k F(z,d/-i) for some matrix valued 
Laurent polynomial L k . 

In consequence, F(z, fd/j,f*) = L(z) + f(z)F(z, dfi)f(z)^ for some 
scalar Laurent polynomial L. This relation proves the proposition. □ 

Although this proposition holds for matrix measures of arbitrary 
dimension d, we will use it only for d = 1,2, which are the cases related 
to QRWs on Z > and Z. Combining this result and Proposition 
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113.11 we see that the local transient states can be characterized by 



\f\ 2 F{z,dfi) e L%(T) in Z > and f{z)F{z,d f i)f{z)^ e L 2 ^{T) in Z. 



These alternative conditions provide a very practical way to de- 
termine the local transient states: they must have as an associated 
function a Laurent polynomial which cancels the singularities of the 
Caratheodory function which are responsible for the non integrability. 

Remember that the state |0) ® |f) of a QRW on Z > is transient 
exactly when F(z,dfi) G L 2 m(T). Then, \f(z)\ 2 F(z, dfi) E L%_(T) for 

2tv 2tt 

any Laurent polynomial /, so we find the following direct consequence 
of the previous results. 

Corollary 13.1. The local states of a QRW on the non-negative inte- 
gers with non trivial distinct coins are all transient if and only if the 
state |0) (8> |t) i g transient. 

Let us apply the previous results to the analysis of the recurrence 
for the examples of QRWs given in Section [TTJ 

Consider first the Hadamard coin on the non-negative integers. The 
corresponding Caratheodory function can be written as 



for some choice of the square root. |F| 2 has two singularities on T: z 
and zq. Neither is Lebesgue integrable, so the state given by a spin up 
at the origin is recurrent. 

Local transient states are characterized by an associated Laurent 
polynomial / such that |/| 4 |F| 2 is Lebesgue integrable on T. There- 
fore, the local transient states are those with an associated Laurent 
polynomial vanishing at z and ~zq. 

Any superposition of up and down states at site has an associated 
function lying in span-fA^X^ = span{l,z -1 } = z -1 span{l, z}. Such 
a function can not cancel both singularities of F, thus any mixed spin 
state at the origin is recurrent. 

However, transient states can appear if we consider a mixing of 
spin states at sites and 1. The Laurent polynomial associated with 
a state a |0> ® 1 1) + b |0) <g> 1 1) + c |1) ® 1 1) is in span{X ,X 1 ,X 2 } = 
span{l, z -1 , z 2 } = z _1 span{l, z, z 2 }, so it can vanish at both, z and 
zo- The Laurent polynomial related to such a state is a + bX\ + cX 2 , 
so it is transient exactly when 




a + &X1O0) + cX 2 (z ) = 0, a + bX^zo) + cX 2 (z ) = 0. 
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Since X\{z) = 1 — \f2z~ 1 and X 2 (z) = —X 1 (l/z) we find that the 
solutions of the above equations are a = and c = —b. This means 
that the transient states with the referred form are spanned by 

|o)®U)-|i)®|T>. 

Following a similar reasoning, and using the form of the third or- 
thonormal Laurent polynomial X$(z) = 1 + V2(z — — \Z2z~ 1 ), it 

is easy to obtain the transient states mixing all the up and down states 
at sites and 1. The result is that such transient states are those lying 
in the span of 

|o)®U)-|i>®|T), |0)®|T) + |i>®U>- 

Let us see what happens if we change the Hadamard coin by another 
equiprobable coin like (I44p . Then, the Caratheodory function 

F{z) = — i^n — 

has a single non integrable singularity at i because — i is a removable 
one. Thus the spin up at the origin is recurrent once again. 

However, in contrast with the Hadamard coin, this QRW has tran- 
sient states at the origin. Such transient states a |0) ® ||) + b |0) ® ||) 
have an associated Laurent polynomial a + bX± which must vanish at 
i. omce Xi{z) = -{y/2z~ l +i) we find that a + bXi(i) = is solved 
by b = i(l + y/2)a, which shows that the states spanned by 

|0)®|T) + i(l + V / 2)|0)(g)||) 

are transient. 

We can also look for the transient states mixing the spin states at 
sites and 1. Using the fact that X 2 (z) = —X\{l/~2) and X 3 (z) = 
1 + iy/2(z — z~ 1 )(i\/2z~ l — 1) we obtain a 3-dimensional subspace of 
transient states a |0) ® |f) + b |0> ® ||) + c |1) ® ||) + d |1) ® \[) given 
by the equation 

a + i(V2 - 1)(6 + c) + (i(V2 - l)fd = 0. 

The QRWs analyzed above are archetypical examples of QRWs on 
the non-negative integers with a non trivial constant coin. They show 
the two possible recurrence behaviours, if one leaves aside the singular 
case of a diagonal coin, which is related to null Verblunsky parameters. 

If the coin is symmetric the QRW is associated with an imaginary 
parameter a = c^ie" 5 . Then, the expression of the Caratheodory func- 
tion F given in Section [9] shows that |F| 2 has two non integrable sin- 
gularities on the unit circle. Indeed, like in the Hadamard case, F 2 is a 
quotient of two coprime polynomials of degree 2 with their roots on T. 
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Hence, the recurrence properties for a (non trivial and non diagonal) 
symmetric coin on the non-negative integers are qualitatively similar 
to those obtained for the Hadamard one. In particular, any state at 
site is recurrent for such a coin. 

On the contrary, a non symmetric coin is related to a parameter 
a = c^ie with a non null real part. The corresponding Caratheodory 
function F has only one non removable singularity Zq on the unit circle. 
More precisely, F(z) = F (z)/(z — z ) with |F | 2 integrable on T and 
Fo(zo) 7^ 0, exactly as for the coin ( 1441) . Therefore, the recurrence prop- 
erties for the coin (1441) are qualitatively the same as for any other (non 
trivial and non diagonal) non symmetric coin on the non-negative in- 
tegers. For instance, these coins always have a 1-dimensional transient 
subspace at site 0. 

Finally, consider the Hadamard coin on the integers. Just as in the 
case of the Hadamard coin on the non-negative integers, the Caratheodory 
function 

has singularities at z = (1 + %)j\pl and Zq, but also at — z and 
— zq. Any of them can cause the non integrability of \fFf^\ 2 for a 
2-dimensional vector valued Laurent polynomial /. 

The local transient states are those whose associated vector Lau- 
rent polynomial / is such that the scalar Laurent polynomial fF^f* 
vanishes at ±z and ±z - On these singularities F is proportional to 
a semidefmite matrix, 

F (±z ) = (l + i)(± 1 F (±z ) = (l-i) 

hence fFof' vanishes on such points if and only if fFo does so. 

Any vector Laurent polynomial / = (a±, a 2 )X + (pi, b 2 )Xi has the 
form f(z) = z _1 p(z) where p is a vector polynomial with degp < 1. 
Therefore, degp-Fo < 3 and fFo can not vanish on four different 
points. This means that any superposition of spin states at sites —1 
and is recurrent. Taking into account the translation invariance of 
the QRW, we find that any superposition of spin states which mixes 
only two contiguous sites is recurrent. 

Thus, a transient state must involve sites which are not contiguous. 
The simplest way to do that is to consider a vector Laurent polynomial 
/ = (ai, d2)X + (b±, b 2 )Xi + (ci, C2)X 2 which corresponds to a state 
P2|-2)(8)||)+6i|-l)(8)|t)+a2|-l)®|J.)+oi|0)®|t)+62|0)®|J.)+ci|l)®|t). 
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Using the expressions 

Xi(z) = (^f 1 , X 2 (z) = X 1 {l/z), 

the conditions f(±z )F (±zo) = f(±z )F (±z ) = become 
(ai,a 2 ) (±i \ + i (t>i,b 2 ) yij +^(ci,c 2 ) = °> 

(ai,a 2 ) f^J +1(61,62) +i(ci,c 2 ) \ ~^^J = 0, 

which have the solutions a\ = a 2 = 0, 61 = c\ and 6 2 = c 2 . That is, the 
transient states obtained are spanned by 

|-2)®||) + |o}®||), |-i)®|t) + |i)®|T). 

Then, the translation invariance permits us to identify as transient 
subspaces all those spanned by states with the form 

I A;) <g> ||) + |fc + 2) <g> ||), |A; + 1)(8) |T) + |fc + 3) ® |T). 

These kinds of results are not specific of the Hadamard QRW, but 
similar recurrence properties hold for any non trivial and non diagonal 
constant coin on the integers. Such recurrence properties are a con- 
sequence of the general expression for the Caratheodory function F 
obtained in Section [101 which shows that F has four singularities on 
the unit circle, with the only exception being the case of null Verblun- 
sky parameters, which corresponds to a diagonal coin. More precisely, 
for an arbitrary non trivial and non diagonal coin, F = (1/ ^/q)F with 
q a scalar polynomial of degree 4 with 4 different roots on T and Fq 
a matrix polynomial of degree 2 which is proportional to a semidefi- 
nite non null matrix on the roots of q. These general results are the 
only ingredients necessary to deduce recurrence properties qualitatively 
similar to those ones obtained above for the Hadamard QRW on the 
integers. 

The fact that any state at a given site is recurrent for an unbiased 
QRW on the integers was proved in [30]. However, the fact that the 
states mixing only two consecutive sites are recurrent too, as well as 
the existence of transient states involving non-contiguous sites is new. 
Moreover, the comments of the previous paragraph show that these 
recurrence properties also hold for any QRW with a non trivial and 
non diagonal constant coin. These general results, together with the 
analysis of the recurrence for QRWs on the non-negative integers, con- 
stitute a novelty which illustrates some of the possibilities of this new 
approach to QRWs. 
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14. Conclusions 

Classical random walks have been traditionally studied using three 
different methods. Two of these have already been used in the quantum 
case. In this paper we propose an approach to the study of QRWs that 
is inspired by the third of these methods. 

Our approach reproduces known results, but also provides new ones 
and new methods of analysis, like the use of the orthonormal Laurent 
polynomials to study the asymptotics or the analysis of quantum re- 
currence using Caratheodory functions. 

This approach can handle non translation invariant QRWs, as well 
as situations where the structure of the one step transitions is richer 
than the ones considered so far. We intend to study some of these 
cases in the future by using CMV matrices where all the Verblunsky 
parameters are allowed to be non-zero; the examples discussed here 
have a simpler structure. This approach can also be adapted in a 
natural way to deal with cases when the walk can go to infinity in 
rather complicated networks as well as in the case of regular networks 
in various dimensions. 

Whereas in the classical case when dealing with an irreducible ran- 
dom walk we have a simple dicotomy: either all states are recurrent or 
all states are transient, we have seen here examples where the situa- 
tion in the quantum case is much more involved. This remains as an 
important area for further investigation. 

It is also important to note that this approach can be easily adapted 
to the case when the number of degrees of freedom in our spins is 
arbitrary. One should consider some of the examples in [30J, such as 
the Grover or the Fourier ones. 

Since the effective use of this approach rests on one's ability to give 
concrete expressions for the orthogonal polynomials and the orthogo- 
nality measure going along with a given CMV matrix there are two 
natural ways to proceed: start with some of the examples where all 
the spectral data is known, such as those in [28], and explore the na- 
ture of the corresponding QRW, or conversely start with some QRW of 
interest and try to compute its associated orthogonality measure and 
polynomials. This is what we have done in this paper. 

It would be nice to look into the QRW that goes along with the 
analog of the Gaussian in the circle, namely the Rogers-Szego case. It 
would also be of interest to study examples where the measure is purely 
discrete. 

These, as well as many other questions, remain as an interesting 
challenge. 



MATRIX VALUED POLYNOMIALS AND QUANTUM RANDOM WALKS 51 



15. Appendix 



Let us calculate the measure fx and the orthonormal Laurent poly- 
nomials x n with Verblunsky parameters a, 0, a, 0, a, 0, . . . for an ar- 
bitrary complex number a with \a\ < 1. For convenience we will omit 
the "hats" in what follows. 

Setting p = a/1 — |a| 2 , the related CMV matrix is 

\ 



C 



a 
P 


p 
-a 










a 
p 


p 
-a 










a 

P o 


p 
-a 















J 

The second recurrence of (J7D, which determines the sequence x,-, 
states that (C — z)x = 0, x = 1, which can be written as 



-z 



a ■ 

,P 
or equivalently 



%2n-l 
%2n 







%2n+l 
x 2n+2 



T 



%2n-l 
%2n 



%2n+l 
%2n+2 

T = -( Z ~ 1 _ ~ a 

p \~a z 



0, 



x = x_i = 1, 



x 



From this identity we see by induction that X2 n (z) = X2 n -i(l/z). 
Using ([9]) and (TTOT) . this is simply a consequence of the vanishing of the 
odd Verblunsky parameters. Thus we only need to calculate X2 n -i- We 
also have 



%2n-l \ _ rpn / I 

Consequently, if \± are the eigenvalues of T, there exist coefficients B± 
independent of n such that 

(49) x 2n - 1 = B+\l + B-\ n _. 
The eigenvalues \± of T are the solutions A of 

(50) \ 2 -p- 1 (z + z~ 1 )X+l = 0. 

Therefore, X2 n ~i(z) are solutions Y n of the second order difference equa- 
tion 

Yn+l + ^n-l = %^n? 



1 



y = 2 ? (* + *- 1 ), 

which is solved by U n (y), U n being the Chebysev polynomials of sec- 
ond kind, given by (j4j). Indeed, the sequences U n (y) and U n -\(y) are 
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independent solutions of this difference equation, thus, there exist co- 
efficients 7, 5 independent of n such that 

X2n-i = 7 U n (y) + 5 C/ n _i(y), 

Evaluating this identity for = 1 and X\(z) = p^ 1 {z~ 1 — a) 

yields 

70) = 1, <*(*) = -p -1 (* + a), 

giving finally 

X2n-i(z) = U n (y) - p' 1 {z + a)U n ^ 1 {y), 1 
(51) y — ( z _i_ 

1 ' X2n{z) = U n {y)-p~\z- 1 + a)U n - 1 (y), 2p l 



This gives the orthonormal Laurent polynomials (xj)°^ . 

To find the corresponding orthogonality measure p we proceed with 
the calculation of the Caratheodory function F(z), \z\ < 1, using ([TBI) . 
Bearing in mind (jUj) we can write 

(52) F(z) = lim X2n -^f> \ z \ < 1, 

"•-* 50 x 2n -i{z) 

where Xj are the orthonormal Laurent polynomials with Verblunsky 
parameters 

—a, 0, —a, 0, —a, 0, . . . 

To take the limit (|52[) we will use the expression (j3HD for Xj, and a 
similar one changing a — ► —a for Xj. The eigenvalues A± of T are 



A± = ^-(z + z' 1 ± V / 0-^ 1 ) 2 + 4|a| 2 ), 

where we choose the square root so that |A+| > |A_| for < \z\ < 1 
(using A + A_ = 1 and A + + A_ = p^ 1 (z + z" 1 ) it is not difficult to see 
that |A+| ^ |A_| for < \z\ < 1 ). Then, (J52J gives 

where _B + is obtained from B + changing a — > —a. 
Using (TJH]) for n = 0, 1 we obtain 
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Thus, the invariance of \± under a — > —a yields 

z~ x + a — p\- z - z" 1 - y/(z - z' 1 ) 2 + 4|a| 2 - 2a 



F z) 



z~ x -a-p\- z - z -i- ^{z - z -i)2 _|_ 4| a |2 + 2a 
y/(z - z- 1 ) 2 + 4|a| 2 + 2Rea _ z - z~ x - 2ilma 



z-z~ 1 + 2ilma y/( z - z ~ 1 ) 2 + 4|a| 2 - 2Rea 

We can obtain the weight w(9) of the decomposition ( 1TTI) for p using 
(I18p by taking the limit r f 1 of Re.F(z), z = re id . Taking into account 
the choice we have made for the square root, when | sin 6*| < |a|, 



lim Jlz - z- 1 ) 2 + 4|a| 2 = ±2\l\a\ 2 - sin 2 9 if J cos ^ - °' 
rTi V V ; 1 1 V 1 1 1 CO s9 < 0, 

while, for I sin#| > lal 



lim J (z - z~ 1 ) 2 + A\a\ 2 = ^2 n /sin 2 9 - lal 2 if J sm ^ - °' 
rTi V V ; V 11 | s i n ^ < 0. 

This allows one to obtain w(6) = lim r |i KeF(re ie ) which is given by 



sin 2 9 



12 



w(9) = -— — - — -, I sin 01 > lal 

v ' \sm9 + Ima 1 ' - ' ' 



and zero otherwise. Equivalently, 



tn\ V S i n2 ^ — Srn2 V n r 1 r 1 

w{9) = . — — , 9 G[rj,7r -rj\U[rj-7r,-rj\, 

| suit/ — sin p\ 

where the angles t] G [0, n/2) and j3 G (—n,n] are defined by 

sin 77 = I a I , sin/? = — Ima, sign(cos/3) = sign(Rea). 

The choice of the sign of cos (3 does not affect for the weight, but will 
be important for the discussion of the mass points. 

Thus the weight is supported on two symmetric arcs of angular 
amplitude 2r\ centered at ±i. 

Concerning the singular part of the measure, its support must lie 
in {e ld : lim^i F(re 10 ) = 00}. From the expression of F{z) given 
previously we see that there is only one possible point in this support: 
e l/3 . Thus it can only be a mass point with a mass given by (jl9l) . which 
yields 

I Rea I yf sin 2 rj — sin 2 (3 

A - \lma\ 2 I cos/?| 



M{e l/3 }) 
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The mass point is located outside the support of the weight because 
| sin/? | < sin rj. Indeed, (3 lies on (—77,77) or in its symmetric arc de- 
pending whether Rea > or Rea < 0. In the limit case Rea = the 
mass point dissapears. 
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